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INTRODUCTION 


A conference  on  Coherent  Radars  for  Range  Instrumentation  took  place  at 
the  Space  and  Missile  Test  Center  (SAMTEC),  Vandenberg  ATB,  California,  from 
S"?  March  1974,  under  the  aegis  of  Range  Commanders  Council  .Exeeutive 
Committee  Member  Mr.  Stan  Radom  (Technical  Director,  SAMTEC).  the  theme 
of  the  conference  concerned  "what  has  been  the  past  performance  and  what  is 
the  future  potential  of  Coherent  Radars  for  test  support  instrumentation."  in 
addition,  emphasis  was  placed  upon  preparations  for  the  upcoming  GEOS-C  launch. 

This  publication  contains  outlines  and  excerpts  from  the  various  conference 
briefings.^  Names  and  addresses  of  conference  presenters  precede  each  of  these 
synopses.  It  is  intended  that  this  document  provide  a means  for  obtaining  addi- 
tional dialogue  and  input  in  the  area  of  Coherent  Radars. 
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SAMTEC  CONFERENCE  ON  COHERENT  RADARS 
FOR  RANGE  INSTRUMENTATION 

5-7  March  1974,  Bldg  7000  Theatre 
Vandenberg  AFB,  California 


0900. 5 March 

Coherent  Radar  as  a Range  and  Range-Rate  Instrumentation  f 'stem  — 

Robert  Green.  White  Sands  Missile  Range  (WSMR) 

Evaluation  of  Range-Rate  Data  - 
Ball  Chin,  WSMR 

Coherent  OuXb  Summary  — 

Bill  Krabill.  NASA 

C-Band  and  TRANET  Tracking  Biases  Relative  to  Collocated  Lasers  — 

John  Berbert,  NASA 

1330. 5 March 

Infiuence  of  Range-Rate  Data  on  ICBM  Instantaneous  Impact  Prediction  Errors  — 
Russell  Roy,  FEC 

Doppler  Track  Evaluation-Operational  Data  - 
Virginia  Fagertin,  FEC 

Coherent  Signal  Processing  Experience  at  SAMTEC  - 
Bill  Collins.  FEC 

Comments  on  SAMTEC  Coherent  Tracking  Data  - 
Maj  Tom  Thomason,  SAMSO 

0900. 6 March 

Probieim  Associated  with  Development  of  a Launch  Head  Range  Safety  System  for 

Containing  ICBM  Launches  in  the  Kwajalein  Laguon  Corridor  — 

Stan  Radom,  SAMTEC 

Accuracy  Requirements  for  Minuteman  Ballistic  Missile  Tracking  — 

Maj  Tom  Thomason,  SAMSO  CLASSIFIED 

Engineering  Improvements  for  Coherent  Radar  Tracking  - 
Renzo  Mitchell,  RCA 


GEOS— C Mission  Objectives/Profile  — 
Ray  Stanley,  NASA 


ix 


1330, 6 March 


GEOS-C  C-Band  Data  Handling  — 

Bill  Krabill,  NASA 

GEOS— C C-Band  Operations/Support  - 
Ben  Jackson,  NASA 

GEOS— C Coherent  C-8and  Transponder  Technical  Characteristics  - 
Alan  Selser,  NASA 

Defe;ise  Mapping  Agency  Test  Objectives  for  GEOS-C  - 
Maj  Larry  Beers,  DMA 

Conference  Summary  and  Group  Discussion  on  ''Where  do  we  go  from  here?" 
0900, 7 March 

Meeting  of  members  of  GEOS-C  C-Band  Working  Group  - 
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GENERAL  INFORMATION 

1.  The  telephone  number  to  be  used  for  recei  /ing  telephone  messages  is  (AUTO VON  276*6190); 
(COMMERCIAL  805*866*6190).  Telephone  i^  ’ ated  in  the  hall  immediately  outside  the 
entrance  to  the  theatre  in  Bldg  700C. 

2.  There  will  be  a military  bus  available  for  and  from  the  Officers  Club  for  lunch. 

3.  Various  telephone  numbers  are  listed  below: 

Scheduled  Airline  Ticket  Office  (SATO),  Rm  C-105,  Bldg  11777:  734-4381 

Vandenberg  Motel,  Lompoc:  736-5605 

Vandenberg  Inn,  Santa  Maria:  922*6631 

Village  Inn,  Vandenberg  Village:  733-3571 

Limousine  Service  to  airport:  736-3636 

Airways  Rent-a-Car,  Lompoc:  736*8521 

Airways  Rent-a-Car,  Santa  Maria:  922*1994 

Heru  Rental  Car,  Santa  Maria:  925-1305 

Army  Liaison  Office:  276*7442/6907 

NASA  Western  Office,  SVAFB:  275*3024 

Base  Taxi:  276*5711 

4.  Mail  can  be  received  in  care  of  the  Technical  Director,  SAMTEC/CA,  Vandenberg  AFB,  CA  93437 
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APPENDIX  A 


COHERENT  rtADAR  AS  A RANGE  AND  RANGE-RATE  INSTRUMENTATION  SYSTEM 


ROBERT  E.  GREEN 
Systems  Management  Division 


WILLIAM  L.  SHEPHERD 
Research  Projects  Office 

Instrumentation  Directorate 

US  Army  White  Sands  Missile  Range,  New  Mexico  88002 


I.  INTRODUCTION 


I 

f 


US  Am\y  White  Sands  Missile  Range  (USAWSMR)  has  nine  coherent  radars  in  its 
instrumentation  inventory.  The  locations  of  these  radars  were  chosen  to 
provide  good  geometry  when  this  equipment  is  used  as  a range  and  range-rate  ; 

instrumentation  system.  The  measurements  of  range  and  range-rate  have  in-  | 

herent  properties  which  are  desirable  for  instrumentation  purposes.  The 
magnitude  of  the  errors  present  in  these  measurements  are  relatively  in- 
dependent of  the  distance  to  the  target.  This  paper  outlines  the  prelimi-  1 

nary  results  of  an  evaluation  of  the  USAWfMR  coherent  radars  as  a range  and  • 

range -rate  instrumentation  system.  Two  methods  of  processing  range  and 
range -rate  date  were  developed  for  this  evaluation.  The  first  method  de-  ! 

veloped  inputs  range  and  range-rate  from  three  or  more  coherent  radars  and  i 

outputs  rectangular  Cartesian  components  position  and  velocity.  The  second 
method  inputs  range  and  range-rate  from  three  or  more  coherent  radars, 
filters  the  range  and  range-rate  data,  performs  a numerical  differentiation 
to  obtain  range  acceleration,  and  outputs  rectangular  Cartesian  components 
of  position,  velocity,  and  acceleration.  This  paper  includes  a description 
of  tne  techniques  developed  and  an  example  of  some  data  processed  using  the 
techniques. 


II.  POINT  ESTIMATION  - THREE  STATIONS 


In  this  section,  the  simple  and  geometrically  pleasing  notation  of  three 
dimensional  vector  analysis  is  used.  The  solutions  obtained  here  are  used 
in  Sections  III  and  IV  to  obtain  start  vectors  for  the  Gauss -Newton  itera- 
tion for  the  N-station  data  processing  techniques. 

Let  0,  A,  B be  three  noncol linear  points  from  which  measurement  of  range  to 
a point  P is  made.  Let  R^,  R2,  R^,  respectively,  denote  the  range  from 

“►  -*■ 

0,  A,  B,  respectively,  to  P,  and  A,  B,  respectively,  the  vectors  from  0 to 

* -*■ 

A,  B,  respectively.  It  is  shown  that  P,  the  vector  from  0 to  P,  is 
P-aA  + bB+cC,  C«Axi  , (2.1) 


♦Larry  A^rml jo,  "Determination  of  Trajectories  Using  Range  Data  From  Three 
NoncoUlnear  Radar  Stations, “ US  Army  Signal  Missile  Support  Agency, 

US  Army  White  Sands  Missile  Range,  NM,  1960. 
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BlAr  ■ g(A.B) 


(R?  - aa  - 

c s ! 

iCi 


a * - r\  + lAl^) 


B * ^ ■'■  1®!^)  • (2.6) 

With  "*  ' 2,  3,  respectively,  denoting  the  coordinate 

vector  of  A,  B,  C,  respectively,  with  respect  to  a^rectangular  Cartesian  co- 
ordinate system  with  the  same  handedness  as  A,  B,  C and  with  origin  at  0, 
(2.1)  can  be  written  as 

+ bap^  + ca^^^  , (2.7) 

’*2  ° **2'^  * *’4^^ 

X3  = aa^^^  + ba^^^  + ca^^^  . (2.9) 


where  x • (x-j,  X2.  Xj)^  1s  the  coordinate  vector  of  P. 

With  pp  p^,  P3,  respectively,  denoting  the  range  measurements  from  0,  A,  B, 
respectively,  to  P the  substitution  of  p^,  pg,  pj,  respectively,  for  , Rg, 
Rj, ^respectively,  lets  one  compute  the  corresponding  rectangular  coordinates 
of  P from  (2.7),  (2.8),  (2.9). 

Differentiating  In  (2.7),  (2.8),  and  (2.9) 
x^  = aaj^^  + ba|^^  + eaj^^ 


(2.10) 


Xg  * , 


Xj  = aa^^  + bal^^  + cBj^^  , 


where  (from  differentiating  In  (2.2),  (2.3),  and  (2.4)) 
4 = i;ii:|2  - 6(a-b)]/|c|2  ^ 

b » [b|A|^  - i(A-i)]/|ci2  . 

c = (2R^R^  - ae  - bB  - ao  - b6)/(2c|Cl)  , 

with  (from  (2.5)  and  (2.5)) 

oi  * “ ^2^2  * 

B - R^R^  - RjRj 


(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 


If  In  addition  to  , P2>  P3  measurements 
available,  (2.10),  (2.11),  and  (2.12)  can 

of  range  rate,  p^,  pg,  P3  are 
be  used  to  compute  x<|,  X2»  X3. 

Differentiating  again. 

•x-  =-a-aSlKb-aS2)+caS3)  , 

(2.18) 

x'2  » aa^'*^  + Ka^^^  + ca^^^  , 

(2.19) 

x’2  “ aa!^*^^  + ta^^^  + , 

(2.20) 

where 


A-4 


(2.21) 


'a  = falBl^  . b(A-B)]/1C|2  . 

b = [6|A|2  -V(A-B)]/|C|2  . (2.22) 

[2{R,’r,  + R?)  - aei  - Zai  - aa  - Bb  - 2bB  - bB  - 2c^lCl^] 

c- LJ ] , (2.23) 

2clCl^ 

and 

a “ R^R,  + r|  - RgRg  - Rj  . (2.24) 

B « R^R^  + R^  - R3R3  - R3  . (2.25) 

With  measurements  p ^ > P2»  V3  also  available ^ (2.18) » (2.19),  and  (2.20)  can 
be  used  to  compute  x.|,  x'3. 

MI.  GAUSS-NEWTON  POINT  ESTIMATION  - N^STATIONS 


In  Section  II,  it  was  convenient  to  use  the  three  dimensional  vector  analysis 
notation.  We  now  change  to  a more  general  notation,  >?ading  to  a systematic 
formulation  of  the  Gauss-Newton  method  for  solving  the  N-station  problem. 

Let  A^,  i = 1.  2 N,  denote  the  points  at  which  the  stations  are 

located  and  AyK  a^^^,  a^^^,  respectively,  the  first,  second,  third, 
respectively,  coordinates  of  A^  with  respect  to  a right-handed  coordinate 
system,  so  that  the  coordinate  vector  of  A.^  is 


Let 


(x^ , X2B  ^3)^ 
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denote  the  coordinate  vector  of  P.  Let 


R 


(Ri 


(3.1) 


where  the  components  of  R are  the  ranges  from  to  P,  given  by 


R,  ■ [ I (x,^  - ■ (3.2) 

(Notice  that  (3.1)  and  (3.2)  give  R explicitly  as  a function  of  x.  We  will 
sometimes  use  R(x)  for  this  function.  Formulas  for  various  derivatives  as<* 
sociated  with  (3.2)  are  given  in  Appendix  A.) 

Three  problems  In  point  estimation  are  posed,  and  their  solution  by  the  Gauss- 
Newton  iteration  described.  A brief  description  of  the  Gauss -Newton  Iteration 
is  given  in  Appendix  B. 

PROBLEM  3.1 

Given  p * (p-|,  P2t  . . . . t a vector  whose  i^*^  component  is  a (noisy) 

measurement  of  the  range  from  A^  to  P>  find  an  estimate  of  x. 

PROBLEM  3.2 

Given  p as  in  Problem  3.1.  and  p.  an  N-vector  whose  i^^  component  is  a noisy 
measurement  of  the  time  derivative  of  the  range  from  Aj  to  P.  find  an  estimate 
of  X and  x. 

PROBLEM  3.3 

Given  p and  p as  in  Problem  3.2  and  p an  N-vector  whose  i^^  component  is  a 
noisy  measurement  of  the  second  time  drivative  of  the  range  from  Aj  to  P. 
find  an  estimate  of  x.  x.  and  x.  ^ 

The  description  in  Appendix  B of  the  Gauss-Newton  estimate  is  phrased  so 
as  to  include  these  three  problems  as  special  cases.  For  each  of  the  prob- 
lems we  first  identify  the  appropriate  vectors  of  the  problem  with  the  vectors 
m.  u.  f.  The  elements  of  the  corresponding  F(u)  are  then  obtained  from  Ap- 
pendix A and  the  start  vector  from  Section  II. 

For  Problem  3.1. 

m ■ p,  u ■ x,  f(u)  » R(x)  . 


A-6 


and  the  start  vector,  Is  obtained  from  (2.7),  , . . , (2.9)  with  R^, 
^2*  R3  replaced  by  pg,  P3  1n  (2.2),  ....  (2.6). 

For  Problem  3.2, 


p i 

< 

R(x) 

. » U ■ 

. • f(u)  - 

R(x,x) 

p i 

and  the  start  vector,  Is  obtained  by  extending  the  above  start  vector 
by  replacing  R^,  Rg,  R^  by  j5g,  In  (2.13),  ....  (2.17)  and  computing 
from  (2.10),  ....  (2.12). 

For  Problem  3.3, 


p 

X 

R(x) 

m = 

p 

, u « 

X 

, f(u)  • 

R(x,  x) 

p 

X 

R(x,  X,  x) 

and  the  start  vector,  Is  obtained  by  extending  the  of  Problem  3.2 
by  means  of  Equations  (2.18),  . . . , (2.25). 

It  1s  worth  noting  that  for  Problem  3.2,  the  objective  function  Involves 

both  p and  p.  Hence,  x,  the  estimate  of  position  for  this  problem  depends 
on  range-rate  measurements  as  well  as  on  range  measurements.  For  Problem 

3.3,  the  objective  function  Involves  p,  p,  and  V.  In  this  case,  x and  x, 
the  estimate  of  position  and  velocity  depends  on  range,  range-rate,  and 
range  acceleration  measurements.  This  attribute  provides  for  additional 
redundancy  In  the  estimation  of  position  in  the  case  of  Problem  3.2,  and 
position  and  velocity  for  Problem  3.3  that  would  not  be  used  If  position, 
velocity,  and  acceleration  were  estimated  sequentially. 


IV.  RANGE  AND  RANGE-RATE  FILTERING  AND  DIFFERENTIATION 


Measurements  of  range  acceleration  are  not  normally  available  from  the 
coherent  radars  at  USAWSMR.  In  order  to  use  the  estimation  technique  de- 
scribed as  Problem  3.3  In  Section  III  of  this  paper.  It  Is  necessary  to 
have  estimates  of  range  acceleration.  This  Information  was  derived  using 
^a  numerical  differentiation  technique.  We  chose  to  dc  the  smoothing  and 
differentiation  on  the  measurements  Instead  of  the  components  of  position 
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and  velocity  to  avoid  correlation  in  the  error  statistics  of  the  data  being 
filtered.  It  also  permits  the  use  of  addUional  redundancy  in  the  point 
estimation  process.  The  filtering  process  developed  makes  use  of  "a  priori" 
information.  In  post  flight  data  processing » it  is  relatively  easy  to  tiiake 
use  of  "a  priori"  information  and  the  improvement  in  data  quality  is  drama- 
tic. In  processing  orbital  information,  the  use  of  "a  priori"  Information 
is  a standard  practice.  Most  of  the  data  generated  at  USAWSMR  is  from 
powered  flight  or  intra-atmospheric  free-f light.  This  makes  the  use  of  dif- 
ferential equations  of  motion  for  "a  priori"  information  very  difficult. 

Since  the  process  we  developed  requires  that  each  set  of  range  and  range- 
rate  measurements  be  filtered,  the  filter  process  is  not  optimal  but  was 
chosen  to  provide  near  optimal  filtering  that  also  provides  computatio  a1 
efficiency.  The  process  developed  during  this  stuc(y  is  a recursive  filter 
which  uses  a cubic  spline  as  the  predictor.  The  "a  priori"  information  about 
the  trajectory  is  Introduced  as  a step  function  for  the  third  derivative. 
Integration  of  this  step  function  is  used  to  derive  a predicted  trajectory. 
This  prediction  is  combined  with  the  measured  range  and  range-rate  data  to 
provide  filtered  estimates  of  range,  range- rate,  and  range  acceleration. 

The  filtering  process  is  described  as  follows. 

Let  Op  J2,  . . . I J|^  be  the  set  of  M steps  of  a step  function  that  approx- 
imate the  third  derivative  profile  of  the  trajectory  under  consideration. 

Let  Tp  Tg,  . . , , Tnn  be  the  set  of  time  intervals  over  which  the  es  - v-  ed 
acceleration  rates  Op  02»  . . . , J|^  are  used,  respectively. 

Let  5t  be  the  time  interval  between  data  samples. 

Let  Aq,  Vq,  and  be  the  initial  values  of  range  acceleration,  range-rate, 
and  range,  respectively.  The  predicted  range  acceleration  is  Jven  b; 

A|^  = + A|^_,  . (4.1) 


The  predicted  range-rate  is  generated  using 

\ + . (4.2) 


The  predicted  range  values  are  obtained  using 


CJ4(6t)-’] 


[Ak.,(6t)2] 


+ V,^.,(6t)  + 


(4.3) 


where  j ■ 1 » 2,  3 M. 
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The  filter  is  mechanized  to  provide  estimates  of  range,  range-rate,  and 
range  acceleration  that  are  a weighted  average  of  the  present  measurement 
and  a prediction  based  on  the  previous  estimate  and  the  acceleration  rate 
profile.  The  values  used  for  the  measured  range  acceleration  are  generated 
using  the  equation 


(4.4) 


The  filtered  estimates  of  range,  range-rate,  and  range  acceleration  are 
generated  using  the  equations 


: . o.(6t)^  - : 


(4.5) 

(4.6) 


J,(6t)^  R.  ,(6t)^  ; 

\ ” ^**3^  * (~^~S ^ Rk.i(«t)  + Rk.i)0  - W3)  . (4.7) 


where  w^,  w^,  are  weighting  values.  These  weighting  values  are  chosen 

empirically  based  on  the  trajectory  being  processed.  The  filtering  process 
Just  described  operates  on  a series  of  range  and  range-rate  measurements  as 
a function  of  time,  while  the  point  estimation  techniques  described  In 
Sections  II  and  III  of  this  paper  operate  on  a collection  of  measurements 
from  several  stations  at  a single  point  in  time. 


V.  RANGE  AND  RANGE-RATE  DATA  PROCESSING  RESULTS 


The  techniques  displayed  In  this  paper  were  applied  to  a set  of  data  col- 
lected at  USAWSMR  on  12  September  1973  (Appendix  C).  The  data  was  col- 
lected by  three  coherent  radars  tracking  a LOKI  Sphere.  This  test  provides 
a very  good  radar  target.  An  aluminized  inylar  balloon  one  meter  In  diameter 
Is  ejected  at  high  altitude.  Radar  is  the  only  instrumentation  available 
at  USAWSMR  that  can  collect  data  on  this  target.  This  limitation  makes 
assessment  of  data  accuracy  very  difficult. 

A segment  of  data  130  seconds  long  was  processed  using  the  range  and  range- 
rate  point  estimation  technique.  This  data  was  encoded  at  20  samples  per 
second  and  then  averaged  to  5 samples  per  second.  No  other  smoothing  was 
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performed  on  this  data.  The  first  20  seconds  of  this  data  contains  some 
errors.  It  should  be  noticed  that  the  errors  are  not  propagated  along 
the  trajectory  using  this  technique.  The  remaining  110  second  segment 
of  data  Is  free  of  ambiguities.  This  110  second  segment  was  processed 
using  the  recursive  filter  to  smooth  the  range  and  range-rate  data  and 
generate  the  required  range  acceleration  data.  A comparison  of  the 
position  and  velocity  data  generated  by  the  two  processes  shows  that 
filtering  of  the  range  and  range-rate  data  did  not  change  these  results 
significantly.  The  acceleration  data  generated  Is  smooth  and  represents 
the  kind  of  accelerations  that  one  might  Intuitively  expect  this  vehicle 
to  undergo.  Visual  comparisons  with  the  range  user's  single  radar  data 
reduction  Indicated  good  agreement.  The  following  table  shows  the 
“a  priori"  Information  that  was  used  to  obtain  the  acceleration  estimates. 


TABLE  5.1  "A  PRIORI"  ESTIMATES  OF  RANGE  ACCELERATION  RATE 


Time  Interval 

Radar  354 
R,  (ft/snc^) 

Radar  123 
Rg  (ft/sec^) 

Radar  352 
R3  (ft/sec^) 

20 

.573 

.720 

.492 

20 

1.746 

2.937 

4.503 

20 

-.990 

-1.701 

-2.178 

20 

-.513 

-.666 

-1.023 

20 

.009 

-.063 

-.138 

10 

.120 

.126 

.018 

Notice  that  the  "a  priori"  Information  used  is  simple.  Evaluations  of 
this  technique  that  have  been  performed  at  USAUSMR  Indicate  that  this 
filter  Is  not  particularly  sensitive  to  errors  In  the  "a  priori"  informa- 
tion. It  appears  that  the  coherent  radars  are  capable  of  producing 
high  quality  trajectory  data  when  they  are  operated  as  a range  and  range- 
rate  Instrumentation  system. 


VI.  SUMMARY 


This  paper  has  presented  data  processing  methods  that  are  applicable  to 
coherent  radar  when  used  as  a range  and  range-rate  Instrumentation 
system.  A technique  for  smoothing  and  differentiating  range  and  range- 
rate  data  has  been  presented.  This  technique  permits  "a  priori"  tra- 
jectory information  to  be  used  easily  and  with  high  computational 
efficiency.  A description  of  point  estimation  techniques  that  can  be 
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used  to  obtain  rectangular  Cartesian  components  of  position,  velocity, 
and  acceleration  from  range  and  range-rate  measurements  Is  Included. 

A description  of  results  obtained  from  application  of  the  techniques 
described  In  this  paper  to  actual  flight  test  data  Is  Included.  Analysis 
of  the  data  obtained  Indicates  that  coherent  radars  are  a potential 
source  of  very  good  quality  trajectory  data  when  used  as  a range  and 
range- rate  Instrumentation  system* 


A-ll 


R -X 

— r-  - \ 

3X. 


bRj 


(A. 8) 


(A. 9) 


3X. 


Differentiating  (A. 3)  with  respect  to  x^,  x.,  and  using  (A. 5),  {A. 6), 
and  (A.9).  ^ 


3Rj  SRj 


2R. 


3R. 
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1 3Xj  3Xj  1 i J 


2R. 


dRa 


3Xa 


•+  R. 


aUa 


3X4 


2X4 


(A. 10) 
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1 axj 


From  (A.1) (A. 12)  successively, 
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(1), 2,1/2 
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(A. 13) 
(A. 14) 
(A. 15) 
(A. 16) 
(A.17) 
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3R4 


3X4 


(A. 19) 
(A.20) 
(A. 21) 
(A. 22) 
(A.23) 
(A. 24) 


From  the  above  formulas,  we  can  form  various  Jacobian  matrices.  For  example, 
with  R written  R(u)  • R(x,  x,  ii),  where  u^  ■ (x^,  iJ , x^)^  and  3R{u)/3x  the 
N,  3 dimensional  matrix  whose  element  In  the  1^**  row  and  column  Is 
3R^/3Xj,  the  elements  of  3R/3x  are  given  by  (A.23). 
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APPENDIX  B 


GAUSS-NEWTON  ITERATION 


A Gauss-Newton  Iteration  general  enough  for  the  problems  of  Section  III  is 
described  without  discussion  of  its  limitation. 

Let  u denote  a p,  1 vector  of  independent  real  variables  and  f,  or  f(u)» 
a vector  valued  function  of  dimension  q,  p < q.  Let  m denote  a q,  1 dimen- 
sional constant  vector.  Consider  the  equation 

m « f{u)  . (B.l) 

u is  a least  squares  solution  if  u minimizes  the  objective  function 

T(u)  • llm  - f(u)||2  » ? (m.  - f.(u))2  . (B.2) 

i-1  ^ ^ 

If  m is  a measurement  vector*  we  will  say  the  u is  a Gauss  estimate  of  u 
for  the  measurement  m. 

The  Gauss -Newton  iteration  is  described  as  follows. 

f(u,  6)  * ||m  - f(u)  - F(u)6l|^  , (B.3) 

where  P(u)  is  the  q>  p dimensional  matrix  such  that  the  element  in  the  i^^ 
row  and  column  is 

3fi(u) 
au'^  3U'* 


With  a given  vector,  let  denote  the  vector  of  least  norm  which 

minimizes  f(u^  6).  We  consider  only  the  case  where 


j(1)  . . 

(B.4) 

whtfe 

(B.5) 
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with  a given  start  vector,  the  Iteration  is 

, 1 « 1,  2,  . , 


If  converges,  seor  to  u,  we  s^y  u Is  a Gauss-Newton  estimate  of  u for 

the  measurement  m.  In  computation,  a way  to  stop  the  iteration  must  be  used. 

With  u^'*^  the  last  computed  Iterate  we  use 
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DESCRIPTION  AND  USERS  MANUAL 
FOR  THE  NR-AM-I  KRTLST  COMPUTER  PROGRAM 


I.  ABSTRACT,  The  KRTLST  computer  program  dumps  and  analyses  radar 
range  rate  data  recorded  at  the  real-time  facility.  It  Is  a very  simple 
program  to  operate  and  requires  no  progranning  experience  to  use  it. 
Anyone  who  has  a legitimate  computer  project  accounting  number  (PAN)  can 
schedule  runs  at  the  WSMR  1108  Computer  B Facility. 

II.  INTRODUCTION.  The  KRTLST  computer  program  was  written  by  0.  P. 
Kroeger  Tl^s trumentati on  Integration  Section)  as  part  of  an  inter- 
organizational  task  effort  by  Analysis  and  Computation  Division* 

Data  Collection  Division  and  Instrumentation  Directorate  to  evaluate 
the  existing  WSMR  radar  range  rate  system. 

The  KRTLST  program  uses  as  inputs  the  data  recorded  In  real-time  at 
the  Real  Time  Systems  Section's  computer  facility  In  Bldg  300.  These 
Inputs  are  data  gathered  by  radars. transmitted  from  the  radar  site  via 
Lenkurt  modems  and  recorded  on  analog  tape  recorded  all  In  real-time. 
Either  In  real-time  (or  In  deferred  time  by  analog  tape  playbacks) 
digital  log  tapes  of  the  mission  data  are  generated  on  the  WSMR  1108 
real-time  computer.  The  KRTLST  program  cannot  be  used  unless  this  real- 
time recording  has  been  range  scheduled  and  accomplished. 

The  KRTLST  program  outputs  setup  cards  for  checks*  radar  data 
samples  of  time*  range,  azimuth*  elevation,  range  rate*  computed  esti- 
mates of  range  rate  from  range  derivatives*  differences  of  range  rate 
and  range  derivatives,  acceleration*  data  one  sigma  error  estimates* 
ard  flags  of  various  predetermined  e**'^ors.  At  the  end  of  the  data  list- 
ing an  automated  written  consolidated  report  of  error  and  reliability 
statistics  is  furnished. 

III.  DESCRIPTION  OF  ALL  STEPS  NECESSARY  TO  ACCOMPLISH  KRTLST  OUTPUTS. 

STEP  NR  1.  Prior  to  mission*  range  schedule  a real  time  re^'nrdlng 
of  all  radars  supporting  the  mission.  Real-time  recordings  are  only  n.ade 
on  missions  which  have  been  scheduled  specifically  as  "real-time  trans- 
mission and  record". 

STEP  NR  2.  Obtain  the  mission  code  (FC  ■ Foxtrot  Charlpy*  AB  ■ 

Alpha  BravoT  et.c.)*  mission  date*  and  mission  start  and  stop  times  If 
available. 
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Obtain  "PSL  Real  Time  Playback  Request"  cards  from  PSl  Input  Office* 


Bldg  300.  Fill  out  request  card  ^see  Fig  (1A))  with: 

1. 

Requestor 

- 

your  name. 

2. 

Acct  Code 

- 

your  accounting  code. 

3. 

Mission 

- 

put  In  mission  code. 

4. 

Program 

- 

put  In  "UNILOG" 

5. 

Date 

- 

put  In  the  mission  firing  date. 

6. 

Phone 

- 

your  duty  phone 

7. 

Building 

- 

put  In  X In  the  300  block. 

8. 

Timing 

- 

put  an  X in  the  range  block. 

9. 

Remarks 

- 

put  "log  all  radars"  plus  mission  start  time  and 
mission  stop  time  If  available. 

10. 

Bottom  Line 

- put  an  X In  the  deliver  box  along  with  your 
name  and  nearest  PSL  delivery  point. 

FIGURE  (1A)  shows  my  request  for  an  1108  computer  log  tape  of  an 
aircraft  track  mission  fired  on  14  Dec  73.  The  mission  was  Charley 
Foxtrot  (FC), 

FIGURE  (1B)  shows  my  request  card  after  the  logging  procedure  was 
accomplished.  PSL  operations  filled  In  the  following: 

1.  Analog  reel  block. 

2.  Tape  classification  block. 

3.  Tape  location  block. 

4.  Tape  reel  number  block. 

5.  Karked  the  ID's  of  the  radars  available  In  the  input  blocks  and 
their  associated  computer  sub*  'Channels. 
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FIGURE  (IB)  shows  that  radar  R113>  R123.  R127,  R352,  R354  and  R373 
were  recorded  and  the  Input  tape  for  the  KRTLST  Is  labeled  W38B/42703. 
With  this  Information  we  are  ready  to  set  up  the  deck  for  a KRTLST 
computer  run. 

STEP  NR  3.  The  KRTLST  run  deck  consists  of  seven  punched  IBM  cards 
and  one  request  card.  (See  FIGURE  (2))  Only  two  cards  have  to  be 
punched  for  each  run.  The  others  remain  constant.  The  two  cards  to  be 
punched  are: 

1.  Card  Nr  1 - Punch  the  Input  data  tape  location  in  columns  16,  17, 
18,  19.  Puncn  the  Input  data  tape  reel  number  In  columns  28»  29»  30,  31, 
32.  For  example,  FIGURE  (3)  shows  data  punch  for  Card  Nr  1 for  the  air- 
craft mission  deslanated  on  the  playback  request  card  shown  In  FIGURE 
(2B),  (W385/42703) 

2.  Card  Nr  6 - The  data  card  Is  the  most  complicated  and  Is  used  to 
control  the  radar  to  be  Jumped,  the  portions  of  the  run  to  be  dumped, 
and  the  ics.9  or  short  print  control  options.  (See  FIGURE  (4)) 

Columns  1 to  4 are  the  easiest.  The>  simply  contain  the  name  of  the 
radar  to  be  dumped.  The  only  allowable  characters  for  the  dump  are: 

Columns  12  3 4 
R 1 1 3 
R 1 2 3 
R 1 2 7 
R 3 5 4 
R 3 5 2 
R 3 9 3 
R 3 5 0 
R 3 5 1 

Columns  11  to  20  contain  the  range  time  of  the  desired  start  time 
of  the  run.  (Three  decimals)  j'Ixample:  67350.000) 

Columns  21  to  30  contain  the  range  time  of  the  desired  stop  time  of 
the  run.  (Three  decimals)  (Example:  68650.000) 
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Columns  59  and  60  are  an  Index  to  be  used  with  Columns  67 > 68 • 69  and 
70  to  control  specific  segments  of  the  trajectory  to  be  dumped  such  as 
first,  second,  third  — n1nety*n1nth  segment. 

Columns  67,  68,  69  and  70  contain  a time  segment  of  seconds  (0000  to 
9999)  to  be  dumped.  (Always  use  in  conjunction  with  Columns  59  and  60. 
For  example.  If  column  60  had  a “3"  In  It  and  Columns  68,  69  and  70  had 
100  In  them,  the  program  would  print  only  the  third  100  seconds  of  data 
on  the  tape.  If  Columns  59  and  60  and  67,  68,  69  and  70  are  used  then 
the  start  time  (Columns  11  to  20)  and  the  stop  time  (Columns  21  to  30) 
should  be  blank.  (And  vice-versa.) 

If  none  of  the  print  time  controls  are  used  t.e.  Columns  11  to  20, 

21  to  20,  59,  60,  67  to  70  are  all  blank,  then  all  data  on  the  log  tape 
will  be  processed. 

Column  60  controls  short  or  long  print.  If  Column  80  Is  "1",  only 
error  flagged  data  and  chock  prints  every  five  seconds  on  good  data  are 
printed.  If  Column  80  Is  "0"  every  processed  sample  (20  samples /second) 
Is  printed. 

EXAMPLES 


Figure  (5A).  Shows  data  card  which  Instructs  the  program  to  dump  all 
data  on  the  Input  tape  for  radar  R113.  (Long  print) 

Figure  (5B).  Shows  data  card  which  Instructs  the  program  to  dump  all 

data  between  times  67350.000  to  68650.000  seconds  of  the  mission  for 

radar  R113.  (Long  print) 

Figure  iSCh  Shows  data  card  which  Instructs  the  program  to  dump 

error  nagged  data  (short  print)  during  the  third  200  seconds  of  data 

recorded  on  the  tape. 

STEP  NR  4.  The  run  request  card  Is  depicted  In  Figure  (6).  As  can 
be  seen,  there  are  three  major  blocks  to  be  f111v«d  out:  Input,  output 
and  Identification.  (Be  sure  and  correctly  label  all  classification 
blocks.  If  they  are  not  labeled,  PSL  will  not  make  the  run.  Do  not  use 
classified  Input  tapes  for  this  dump  program  because  the  outputs  or  the^ 
oT the  program  will  oecome  classified.  As  a rule  run  only  data  from 
unclassified  projects. 


In  the  Input  tape  block  two  tapes  have  to  be  entered.  The  program 
tape  Is  always  E725/39027  until  after  Feb  74.  Call  Kroeger  (678-1620) 


for  new  program  tape  designations  and  at  four  month  Intervals  after  that. 
The  second  tape  will  be  the  Input  data  log  tape  you  have  requested  from 
the  PSL  Real  Time  Playback  Facility.  The  number  of  cards  will  always  be 
seven. 

In  the  output  section  (under  PRINTS)  always  put  N0RM  unless  you  have 
changed  the  deck  to  get  one  or  more  carbons. 

In  the  Identification  block  check  R under  option  block»f111  In  w>ur 
RUN  IDf  PAN»  requestor's  name*  phone*  organization*  check  Computer  b* 
put  your  delivery  destination  point.  The  time  estimate  and  page  estl* 
mate  Is  computed  as  follows: 

TIMING  ESTIMATE 


1.  Four  minutes  for  first  0 to  200  seconds  of  flight. 

2.  Add  1.1  minutes  for  each  additional  50  seconds  of  flight. 

3.  Round  fractional  times  up  when  equal  to  .5  seconds  and  greater. 
EXAMPLE  NR  1 - For  a 556  second  flight. 

556  sec 

200  - first  seg 

52J55?  7 * 1.1  - 7.7  - 8 


TOTAL  TIME  « 12  MIN 


I 4 min 


EXAMPLE  NR  2 - For  a 320  second  flight. 


320 

200 

1 4 min 

TZff 

2 

SOjTZff 

2 X 1.1  ■ 2.2  . 2 min 

TOTAL  TIME  ■ 6 MIN 


1 


i 


EXAMPLE  NR  3 - For  a 1298 
1298 


TIME  » 4 + 


PAGES  OF  PRINT  ESTIMATE 


second  flight. 

I 4 min 

21  * 1.1  « 23,2  • 23 


23  « 27  MIN 


Ten  pages  plus  one  page  for  every  2.5  seconds  of  flight  data. 
EXAMPLE  - For  1300  seconds  flight. 


10 

520 


?35 


pages 


520 

^.5lT3B8g 


530  PAGES 

IV.  INTERPRETATION  OF  OUTPUT  LISTING.  The  OiJtput  listing  has  four  parts. 
The  frr^  part  Is  a three  page  computer  summary  of  loading  parameters 
which  Includes  the  tape  assignments  for  checks.  The  second  part  Is  a 
one  page  listing  printed  by  the  program  giving  the  data  card  setup  and 
the  error  code  definitions.  The  third  part  (the  longest)  Is  the  listing 
of  the  data  as  follows  (left  to  right) 

1.  Radar  name  or  the  error  code  If  an  error  has  been  detected.  The 
following  list  gives  error  code  definitions. 

ElOOOOO)  - DOPPLER  DVES  VALID  FLAG  HAS  CHANGED. 

E 20000)  - DOPPLER  21PT  RAW  DATA  SIGMA  HAS  CHANGED. 

E 3000)  “ VELOCITY  EEOR  G.T.  14.31  F.P.S. 

E 400)  • TIME  HAS  A GAP  (NOT  .05  SEC)DERIVATIVES  BAD. 

E 50)  • RADAR  HAS  CHANGED  ITS  TRACK  MODE. 

E 6)  - NOT  PRESENTLY  ASSIGNED. 

E 20050)  • ERRORS  2 AND  5 BOTH  PRESENT  (EXAMPLE). 


\ 
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2.  Time  - seconds. 

3.  Time  difference  • seconds  for  missing  data  checks.  (This  may  be 
changed  later  to  display  AGC. 

4.  Radar  track  mode: 

S • skin  track 
B * beacon  track 
N • no  track 

5.  OVES  valid  flag: 

1 ■ valid  flag  yes 

0 « valid  flag  no 

6.  Dopple**  skin  return  flag: 

1 ■ yes 

0 ■ no 

7.  Doppler  COHO  beacon  flag: 

1 ■ yes 
0 • no 

8.  Range  - feet. 

9.  Azimuth  - degrees. 

10.  Elevation  - degrees. 

2 

11.  Second  derivative  of  range  (acceleration)  51  pt  filter  ft/sec  , 

2 

12.  Doppler  velocity  derivative  (acceleration  51  pt  filter  ft/ sec  . 

13.  Range  first  derivative  (velocity)  one  sigma  error  estimate 
ft/ sec.  (21  pt  filter) 

14.  Raw  doppler  velocity  error  estimate  (21  pt  filter)  ft/sec. 

15.  First  derivative  of  range  (velocity)  21  pt  filter  ft/sec. 

16.  Raw  doppler  velocity  ft/sec. 


17.  Velocity  error  (difference  between  raw  doppler  velocity  and  51 
pt  range  derivative)  ft/sec. 

18.  Velocity  error  (difference  between  raw  doppler  velocity  and  21 
pt  range  derivative)  ft/sec. 

The  long  print  (data  card  COL  80-0)  gives  the  above  information  for 
each  and  every  point.  (20  SAM. SEC.) 

The  short  print  (data  card  COL  80»1)  prints  only  flagged  errors  (see 
1.  above)  plus  24  consecutive  samples  prior  to  a flagged  error  and  24 
consecutive  samples  past  the  same' flagged  error.  Also,  during  long 
stretches  of  good  data,  a checkpoint  will  be  printed  every  five  seconds. 

The  fourth  part  of  the  listing  is  a written  statistical  report  of  the 
entire  run.  The  basic  criteria  for  most  percentages  and  reliability 
estimates  are  based  on  the  premise  that  when  the  radar  is  in  track  (either 
skin  or  beacon),  there  should  be  valid  unambiguous  velocity  data  (with 
valid  vel  flags).  The  run  printout  will  not  commence  until  the  radar 
has  been  in  solid  track  for  at  least  two  seconds.  The  radar  reliability 
number  is  based  on  outages  versus  skin  or  beacon  track  once  the  run  has 
begun. 

If  the  computer  run  aborts,  check  the  last  page  of  your  listing  for 
such  things  as  tape  problems,  max  time,  max  pages.  If  either  max  time 
or  pages  Is  indicated,  your  estimates  on  your  run  request  card  are  wrong 
and  should  be  Increased. 

V.  PROGRAM  LISTING.  i*or  those  who  wish  to  know  what  the  program  KRTLST 
coding  Too  ks  nice",  a listing  of  the  program  is  included. 

For  those  who  wish  to  have  deck  of  the  program  for  their  own  purposes 
may  easily  do  so.  The  Fortran  elements  are  on  the  program  input  tape. 

To  punch  (on  line)  the  deck  during  a run  place  the  following  cards 
directly  after  the  copin  card  in  the  data  deck: 


1. 

8PCH,SC 

TPF$.R0YSA 

2. 

@PCH,SC 

TPFS.ROYSO 

3. 

8PCH,SC 

TPFS.LOG 

t 

+ 

COL  1 

COL  14 
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REAL  TIME  PLAYBACK  REQUEST  CARDS 


SAMPLE  DECK  AND  REQUEST  CARD 


FIGURE  (2) 
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INPUT  TAPE  AilGN  CARD 
FIGURE  (3) 
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THE  DATA  CARD 
FIGURE  (4) 


DATA  CARP  (PROCESS  ALL  DATA  - LONG  PRINT) 
FIGURE  (5A) 
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>u> 

'2t 

A3» 


16NHE«T 


t<pe»  relative  location,  NAHtI 


‘GVO  lOL 
1327  |2Sk 
m<t*  30L 

*017  I A 
EiN  1 1 
001 

>Q7<4  tVALlB 
2»B  L 
2%0  HORAO 


OOqI 

000  t 

0001 


0003*3  noL 

OOOol^  I&l** 
0000^3  3^V 


OOOH  \ OOOOO*  UOh-‘B 
0000  t 00*2*3  IJBK 
0003  1 000133  I«0 
OQOS  I OOOODO  IVALIO 

0000  X 00*23H  LCOOl 
QOQN  I OOOoOI  PARITY 


OOOt  00027* 
0003  00030T 
0000  00*2** 
0003  I 000171 
OOOH  I 000003 
0000  3 00*1*3 
OOOO  \ 00003* 
0000  \ 00*2** 
000*  R OOOOQl 


, I*L 

mt 

LD 

iLlft 

\ftSTAt 

IHAVU 

h 

SRl> 


ODOl  OOQl^*'  HOL 

QOQl  OOOOl*  *0L 

' OOOM'T  000000*^071 
0000  I 00*2*3  lOT 
OOOO  00*^77  INjPt 
0003  1 0002*S  IITAT 
QQOQ  r 00*267  I*C 
0000  1 00*2*3  MASK 
■ QOOi"*^  OOlJOOO  T 


SURROUT  [NE  t'NPfcCK 


CUMHON/NRlC/KK,SKtP 


OlHENStON  fAPRECl32*»n 
integer  ENDPIL, parity 
data  N/O/PARITY/0/ 
data  ioNES/O777777aOO0OO/ 
data  Ha*R/07777H0000000/ 
17  IER071L»nE,0>  60  TO  10 
ir  |AAp»ITY«NC,01  60  TO  I* 
ir  <N*gT,0)  6o  to  3* 

10  ENOriL  • 0 

call  NTRAN  I2,I0I 

|6  PABITT  ■ 0 


ui*  data  eiTS,  One  51«N 
OCT  13 


^POSITION  TAPE  TO  LOAD  POINT 


Vll 


*1 


APPENDIX  6 


CALL  NtRANJ2*2»32**1|TaPRECiLCCDE) 
CALL  NTKAN(2,22» 

IK  (LCOOE.GT.OI  GO  TO  30 

IF  (LC00EtEQ.-3)  GO  TO  1 10 

IF  (LCoOEtEUf-2)  GO  TO  120 

30  ICONTR  ■ BOOL  I TAPREC< I ) * 

L»2 
M a 0 

35  N ■ N «.  I 

IF  (N»gTi20)  GO  To  20 
P :»FLOt 18  I 18  ,TaPR£C<L  ) I 
ILIFT*  FL0l0i6iTAPKEC(Ln 
NORAD  • t IWC«1  >/S 
T«TAPREC(L*1 )/l000# 

DO  60  !■!  fNORAO 
IRSTAT  • BOOL < T APREC I I ) 
lOT  ■ rLDI 20.7  » IRST*T  ) 

II  ■ I0T*l 

IF ( 1 I .GTt30>  I I • 30 
IM  I0T.EQ**iH)  I I«M 
IF(  IDT.EQ.MOI  II»Z5 
1F( lOT.EQ.Hl > I 1*26 
IF i lOT .EQ.M2 ) I 1*15 
IF  I 1DT*£Q«<47  > I 1»16 


OREAD  1 Physical  record 
brecord  ok 

UTCST  FOR  bad  record 
8TEST  FOR  EOF 

©control  0»RA0 » 1 ■met  .2«RAd^HET 


©NO*  OF  RADARS  PER  logical  R(-C 


U RADAR  39^ 

© R3S0 
« R3B1 

0 radar  3S2 

(jRADAR  393A 


C 

C ISTAt  is  a new  status  WORD*  example  02S007010330,  0250  • 

C AdT*SYNC,  Q7  • CONF,  counter.  01  ■ run  no*.  03  • T6T  nO». 

C 30  ■ SUBCHANNEL  NO* 

c 


FLDII2.8.IRSTAT) 
FlD(6 .6, I RSTAT } 
FLD(0»3. IRSTAT) 
FLOU  .3,  IRSTAT  ) 
Fl0(6. 6, IRSTAT) 


FLDI R »e* ISTAT (ID) 

FLDI  12.6.  ISTATI  ID) 

FLDI2I  .3.  ISTATI  ID) 

FLD(27|3.  ISTATl  ID) 

FLOI  30,6.  ISTATU  1 ) ) 

I0( I D ■ lOT 

IRIID  • FLDD2 ,2H  .TAPREC  IM  + 5)  ) 
lADD  ■ FLO(  12, 17  ,tAPRECIM*6)  ) 
lEDD  ■ FLDD2.17,TAPREC1M*7)  ) 
IRO(II)  ■ FLD(  13.23, TAPREC(M*0D 
ITDDD*  FLD(33,3.IRSTAT)*16 

ivalioi I D ■ 0 

I VAlIDD  I )"FLD<  32,  I .IRSTAT) 

I VALIS ( I 1 ) "FLOl 31 , 1 , IRSTAT ) 

1 VAlIBU  D«FLDI  30. 1 , IRSTAT) 

1F(  I I *GT*8)  GO  TO  1)8 
60  TO  121 
1 1 B CONTI NUC 

lF(FLD(12il,rAPREClH<>8D*EU*D 
BR  F0RMAT(2( 1X.012) ) 


© RADAR  ID 
Q RANGE 
e azimuth 

9 ELEVATION 
© RANGE  RATE 
©N0V73  FOR  MPS36 
0 N0V73  FOR  MPS36 

« DOPPLER  SKIN  flag 
8 DOPPLER  BEACoN  FuAG 
M MPS  36  NOV  73 


1RD(  I D«  1R0(  I D • i »D 


Il9  continue 
GO  TO  128 


C FPS  16  N£G  ROOT  NOV  73 

121  continue 

1F(FLDD2.1.TAPREC(H«8)))127,)28.127 
127  !RD(  I !)■  IROI  I D'M 

I Ro  D D * OR  I mask  . 1 ro  n D ) 


B-23 


FP5  16  NE<» 

120  COMINUC 

TIMt  SKIP 

|F(TtUT*5KIrl  (iO  To  3 

I 10K«  FLOIO,3a  , If^STfcT  ) 
lFm.tW.KKI  aRITEI6,0H)  llBKi 
3 continue 

H ■ 

60  continue 

M ■ H*2 
80  L ■ L*I«C*1 
RETURN 

1 10  .'ARITY  ■ LCOOt 
call  NtRAN  < ? i22 I 
RETURN 

120  ENqFIL  • I 
RETURN 
END 


HOOT  NOV  73 

« mOUSlKElPI^Q 
I I0K  housekeeping 

B SET  flag  FqH  PakITY  eHROW 
W SET  end  of  file  flag 


COMPILAI lON! 


NO  oiaqnostks 


)ONU 

IKAO 

ooo** 

I 

000133 

IRD 

OOOO 

1 

001132 

IRSIG 

oooo 

I 

1 RUA 

ISH 

IRUN 

oooo 

I 

001  13S 

ISPTS 

DOOR 

00026S 

istat 

000^ 

-f- 

■UB0Ul2 

1 TO 

202 

itest 

0003 

1 

00321N 

ITIME 

doom 

1 

000227 

ITQ 

I 

00001 1 

ITR 

217 

lUM  10 

0003 

I 

00321S 

IUM6 

0003 

1 

003216 



0006 

I 

00007N 

1 VAU  1 

i03t 

IVALIS 

0000 

I 

001  130 

I VNEG 

OOOO 

1 

TvPos 

oooo 

I 

001166 

J 

1S7 

K 

OOOD 

1 

00)160 

L 

oooo 

>■ 

-<»tni63 

nPoint 

0006 

I 

OOOOOl 

PARIT 

1B2 

RTOG 

0000 

R 

001  ISI 

scale 

R 

001212 

SOREL 

0003 

000032 

SDVL< 

>027 

SIGR2| 

0000 

R 

001 167 

oooo 

R 

000223 

SPRI 

OOOO 

H 

001213 

SRREl 

176 

T 

0000 

R 

001JJ>- 

-Test 

oooo 

R 

001207 

tf 

OOOO 

K 

001162 

TI 

lONH 

tlift 

0003 

> 

«^»tyffo3s 

tstart 

0003 

R 

000Q36 

TSTOP 

0003 

K 

000021 

VEL 

. 2Q6 

vlpr 

R 

001216 

VPP 

OOOO 

R 

001216 

XPER 

OOOO 

R 

001166 

XRAW 

)000 

COMWON/MaIn/  IER»  IEt4»i  TIN,DT,ITR,ITDJD>.IT03iRiAiE|VCLi 
tORBl iDk21 •VEL21 iODR&I iDVEL^I »S1GR21 f i , S VEL2  I •S0VL21  » 
2DEL5liOEL2l|T5TART,TSTOPiINAM,lRAi),|^':;'«E*3>|TLirT,OELTi 

3fiurr  ( 20  |6  1 ) I 1 NAME  (9  ) I I chan  ( 9)  , !▼  t M:.  , JUMAi  1 uH9  i tUHiO  » lOPT 
COMMON/ INT/Z I I R ( 30) I 1 A ( 30) I IL (30 } I iNM(30) • 10 ( 30 ) I |TQ (30 ) t ISTAT (30  ) 
COMMON/EOFF/ENOriL#P*RlTY,  |COfJTR|  ILIFT 

COHMON/FeRATE/tVAL 10(30) i 1 V AL I S I 3U ) i 1 V AU B ( 30  ) 

DIMENSION  FINVARIft) ) 

DIMENSION  HOlD(2i‘4|B)  I 1PT(2) 

DIMENSION  SPRKlBt?^) 

INTEGER  ENDFILiPARiTY 
data  iMODC/lHBiIHSilHN/ 

DATA  InAME/RHRI  13fMHRl23iNHRl2  7,MHH39N|qHR3S2t<4HR393lHHR3fi0|NHK35i 
I iNH  / 

data  1CHAn/2i5i8 t in  I 1 5 . 1 6 1 2B » 26  * 30/ 
lPOSPT-0 
INEfiPT"0 
IVPOS  *0 
IVNEG  >0 
I0S1G«0 
1RSI6-0 
INNN  bO 
|NPTS*0 
ISPTSbq 

IBPTSbg 

IDBF 1-0 
1DQF0«0 
1KP*0 
IPTCEbq 
IPTC«0 
IPTCBbO 
lOOPlBO 
|O0P2*0 
IPTU  )*2l 
IPT ( 2)bS1 

RE AO (5, 2002)  I N AM  , | UM6 , TST ART , TSTOP , CELT i Tu I F T 1 1 U« 1 0 i IT1MEiIUm9i 
I lOPT 

WRITE  (6i213) 

WRITE  (6|21N) 

WRITE  (6,212) 

WRITE  |6,2IS)  INAM 


a COUNTER  OF  POS  VEL  PoINT3(Gt  |00  FPS) 
0 COUNTER  OF  NE6  VEL  PQlNTSILT  lOO  FPS) 
e VALID  POS  DOPP 
« valid  NEG  DOPP 


u raqar  intkack  (must  be  gT  So  for  program  To  run) 
w numder  of  points  processed 
0 counter  for  Doppler  sicru  reT'urn  points 

U counter  for  DQPPLER  COHO  BEACON  R£TuRN  POINTS 

9 data  Is  bad  BuT  Fug  says  good 

9 DATA  IS  GOOD  BUT  FLG  SAYS  BAD 


B-2S 


WRITE  TSTART 

WRITE  TSTOP 

WRITE  oeiT 

WRITE  16,212) 

WRlTE(6i220) 

wRITEI6i22n  lUHlOilTlME 
WR1TEI6|212) 

IFU0RT*EQ»0)  write  16,223) 
in  I CRT  • EH  • I ) WRITE  (6,222) 

IF ( IORT*LT*OiOR* I 0PT#6T* I ) WRlTEI6*22*n 

3000  fORHAT(S5H  The  following  list  OIVES  error  code  OE^|N|TIONSt 

3001  rORHATlSSM  £(00000)  ■OOPFUER  OVES  VALID  FU*G  HAS  CHANGED* 

3002  rORMAT(55H  E 20000)  -POPPLER  21PT  RAW  DATA  SUMA  HAS  CHANOEq. 

3003  rORMATlSSH  E 3000)  ■VELOCITY  ERRqR  6*T*  l*),3l  F«P*S* 

300*)  F0RMAT(55H  E MOO)  ■TIME  MAS  A GaPINOT  ,05  SEC  ) DER  1 V AT  1 VES  BAD 

3005  F0RMAT(55H  E 50)  -RADAR  HAS  CHANGED  I TS  TRACK  MODE* 

3006  FORMATlSSH  E 6)  «NOT  PRESENTLY  ASSIGNED* 


3002  r0RMAT(55H  E 20000)  -POPPLER  21PT  RAW  DATA  SUMA  H*5  tH«N«*.g* 

3003  rORMATlSSH  E 3000)  -VELOCITY  ERRqR  6*T*  ^‘♦•3l  F*P*S* 

30QR  F0RMAT(55H  E MOO)  -TIME  MAS  A GaPINOT  ,05  SEC ) DER 1 V AT  I VES  BAD 

3005  F0RMAT(55H  E 50)  -RADAR  HAS  CHANGED  I tS  TRACK  MODE* 

3006  FORMATlSSH  E 6)  -NOT  PRESENTLY  ASSIGNED* 

3002  FORMATISSH  ^ ^ 

3000  F0RMAT(5BH  E 20050)  -ERRORS  2 AND  5 BOTH  PRESENTIEXAMPlC ) 

WRITE(6|2I2) 

WRITE(6«3000) 

wRIT£(6i3007) 

WRITE(6,300I ) 

WRITE(6|3002) 

WRlTE(6i3003) 

WRITE(6*300R) 

WRITE(6,3D05) 

WR1TE(6|3006) 

WRITE(6*3007) 

WR1TE(6i300B) 

WRITE(6i3007) 

WR1TE(6,2I3) 

WRITE  (6i225) 

WRITE  (6i226) 

WRITE  (6,227) 

WRITE  (6,3050) 

WRITE  (6,225) 

C FROM  data  CAROS  COMPUTE  RADAR  ID  FOR  UNPACK  AND  SET  FL AOS • CONST AnTS 


DO  1 t-1,9 
11RR*1CHAN( 1 ) 

tF(lNAHE(|)tEQ*|NAH)0O  TO  3 
CONTINUE 

FORMATIffH  THE  RADAR  SPECIFIED  ON  THE  DATA  CARD  Is  NOT  CATALOGED  a 
IS  A RAN5E  RATE  RADAR  IN  TH 1 S PROOffAR,  STTH*'* ) 

WRITE(6,212) 

WR|TE(6,212) 

WRITE(6,2) 

WR1TE(6,212) 

WRIT£(6i212) 

STOP 

CONTINUE  _ ^ 

SCALE  ■ 3*IRI592653S09793OO/I2*9«V»651 V-ANCtr-nr^RADTANS 

WTOG  ■ 67, 2957795131  B ANGLE  IN  DEGREES 

lRUA-TSTART-100O**TST0P*tO0a* 

IRUN-IRUA  ♦ lUMlO^ITlMC 


-W6 


lRUB«lyMlO*  ITIME 


c Shift  data  holding  arrays 
c 

100  DO  ZOlttltB 
DO  10J«ltRQ 
K«  J*| 

10  0UFF1 : I J) “BuFr  I I I K » 

20  CONTINUE 

00  N0I*|,9 
K«i  ♦e 

DO  >0J«1|N0 
L*J*I 

30  Burr (K  I JlaBuFF (K |L) 

NO  continue 


C LOAD  INRUT  values  FROM  REAlTIML  LOG 
c 

call  unpack  w buffi 

c 

BUFFIS.ei >"  -too-o  0 

ITQI  I IRR)>1TQ(  t IRRI-U 
1F(  ITQUIRR)  tCQiN  IBuFFIBtBl  )*HiO 
IFIITQ ( 1 IRR ) vEOtS  IBUFF IS |6 1 laStO 
BUFFI  1 ,81 l»tR< 1 IRR) 

BUFFIl  ,61 )*BUFFI 1 fSl  )*3.0 
BUFF  12 ,81  I >|RD( I IRR I 
BUFF  I 2, 6 I >«BUFFI2f81 )*i007S 
BUFFI3,ei >•! A I 1 IRR) 

BUFF  I 3 |8t >«BUFF I 3 *61 ) • SCALE 
BUFFIN  ,81 1»1LI IIRR) 

BUFF{H,ei laBUFFINtei }•  SCALE 
Z"Z-*090 

BUFF ( 1S,N|  }«Z-2«0 

BUFF(8,61 )■! VALIDI I iRR 1 

BUFF  17 ,81  )«BOOLU  VALtSi  1 IRK)  } 

Buff  I a 101  i«booli iv alibi i irr)  i 
c 
c 

IF(  irun*eq*o)  go  To  nho  e 

iFi  irub«eq*o)  go  To  4mi  e 

IFIIUM10.LT.2I  GO  TO  NN^  9 

ITST  •IlUM10-n»(lTlME*i:0) 
tOOPl-IOORUl  0 

IFII00P1.LT.1T5T)  GO  TO  100 
1RUN«0 
60  TO  NHO 

HNl  IF  I8UFFI |5,N I ) .LT.TSTART)  GO  TO 


i8l  } > NLWEST  POjNT 

track  mode  negative^  no  TRACK 

0 TRACK  mode  is  beacon 

0 track  mode  is  skin 

0 range  yds 

0 RANGE  FT 

6 range  RATE 

0 range  rate  in  F<PtS* 

0 AZIMUTH 

e azimuth  in  radians' 

0 ELEVATION 

Q elevation  jn  Radians 


IF  transfer,  process  whole  TAPE 
IF  TRANSFERiUSE  TsTART  8 TSTOp 
IF  TRANSFER. PROCESS  1ST  ITImE  pTS 

space  Tape  over  ist  n itihes 


100 


IRUNaO 

HNO  continue 

c 

c the  following  test  does  not  permit  run  to  begin  until  data  arrays 
c HAVE  Enough  in*track  radar  data  to  initialize  Fjlters 

C ISl  consecutive  points  in  beacon  or  skin  ThACK) 

c 
c 

IFI INNN«GT.S0)GQ  to  1002 

B-27 


c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

iTlC« 

I 

(. 

t 

I 

I 


INNN  »o 

DO  1001  1*1 |5t 
I^15 

irtBurriBij) faTfO.oi  innn«innn^i 

1001  continue 

Tl*  BUrFIl&iNI)  Q START  TIhE  OF  FUN 
GO  TO  100 

1002  continue 

SMOOTH  range  and  range  rate  using  21  AND  51  POInT  SnOOThINq 

00  HI  :-li2  » fit  2lPT  THEN  51PT 

NPOI NT  • IPT  n » 

DERkEY  •1*0 

DO  H2  J»li2  e FIT  RANGE  THEB  H •GEKATE 

DO  H3  k"1iB1  « LOAD  RAW  DATA  (61  PT  ARRAY) 

H3  FINVAR(K»"BUFF( JiK) 

call  LLL  LLS  is  an  unconstrained  least  SQUaRESi  2nd  ORpER 
MOVING  ARC  filter*  IT  IS  COMPLETELY  COMPaTaBLE 
THE  The  standard  wsmr  data  reduction  vaa  program* 

CALL  LLL(NP0INT|DeRKCY|FtNVAR,XRAW»xSM,5IGX*0XSM»QS|GX,D0XSM| 
IDDSlGXiCVAR) 

HOLD ( I lU 1 1 ) "XKAtt  b 

H0LDUtUi2>«XSM  W 

HOLD! 1 ,Ui3>*DXSM  b 

HOLOU|J|H)«DpXSM  e 

HOLD! 1 ,JiS)bSQKT(SiGX)  b 

H0LD( 1 iJ«A>bSORT(DSIGX)  9 

HOLD( I tUi7)«50RT(DDStax)  b 

HDLD( 1 ,UtB)«CVAR  9 

H2  CONTINUE 
HI  CONTINUE 

LOAD  INTO  OATx  buffer  ALL  SMOOTH  DATA,  FlAGSi  AND  VARIANCES 

THAT  are  To  be  saved  For  a two  second  interval  •••used  for  print/eoit 


Raw  mid  point  value 
Smooth  mio  pOint  value 
SMOOTH  MID  point  1ST  DERIVATIVE 
Smooth  mid  point  2nd  derivative 
smooth  hid  point  1 SIqMA  err*  est* 

1ST  DERIVATIVE  i SIcHA  ERR*  EST, 
2nd  DERIVATIVE  1 SIGMA  ERR*  EST, 
Variance  or  thE  total  data 


BUFF (9,H1 )•  H0LD(2i 1 »H) 

BUFF ( 10,H1 )kHQLD( I ,2i3) 

BUFF( 1 1 ,H1 )■HOLD( 1 i 1 i6> 

BUFFI 12|H1 )«HOLO(  1 ,2 *51/ *3279 
BUFFI 13|H| )aH0L0(2,l  ,3) 

BUFFI 1H,H1 )•  mold! 1 |2i2l 

test  ERROR  PRINT  OUTS  AND 


Q RANGE  2ND  DERIVATIVE  BIPT 

0 Range  Rate  derivative  2IPT 

9 RANGE  1ST  OERIV,  SIGMA  21PT 

e Range  rate  sigma  2ipt 

0 RANGE  1ST  OERIV,  SIPT 

9 RANGE  Rate  (smooth)  2ipt 

SET  flags 


lEin-o 

THE  TEST  for  Equality  between  NON<*lNTeGERS  may  not  be  meaningful, 
IFIBUFFI A,HO) •NE*BUFFI6*H1 ) I lEl M«100000  H HAS  OOP  FLAG  CHANGED 
IEI2)«0 

tF(BUFFU2,Hl)  ,GT*  1 *2)  IEI2)  *20000  BvEL  SIGMA  6T  1,2  FPS 
!EI3)«0 

TE$T«ABS(N0L0(2i 1 •3l*H0L0( 1*2,1)) 

|F(TEst*«T,  IH,31  ) ICU)*3000  (BVEL  ERROR  GT'IH*3I  FPS 

ICIRI-D 

B-28 


vm 


TEST*  AbS<bUFF  ( l&iH  I )-BUf  F n & I'^O)  ) 

!r(TC5Ti6T,,06,OR«TtST#tT..O‘*bl  |fe(t<)«HOO  it)  TIME  6aP  tRHQW  TEST 

te(S)>o 

bT{C*  The  test  FOR  EQUalMY  HETAEEN  NOn-INTEoERS  may  NOT  Be  MEANlNOFULt 

IF  (BuFF (StHC ) tNE«HuFF (SiHl M lE<S)*SO  w H*0AM  ChAnbEO  TRK  hU 

NOT  useo 
c 
c 

c flags  for  crrohs  frinteo  on  The  Output  listing  in  the  fjRst 
C column  ANO  mill  AFFCAR  as  follows  <0R  ]N  COMBINATIONS): 


IERR«Q 

IERR»  Um*lE(2)4>lEl3)*lE(«t).»IE(SI«IE(6) 
BUFF) 17»MI I • {ERR 


c 

LIOOOOO  ) 

■ 

OOPPLER 

TRACK  Flag  has  Changed 

c 

E020000  ) 

m 

DOPPLER 

1 SIGMA 

G*Tt  1*2  FPS 

c 

1003000  ) 

m 

VELOC I TY 

ERRqR  6 

*T,  1H.3I  FPS 

c 

EOOCROO  ) 

m 

time  has 

A GAP  - 

DERIVATIVES  BAD 

c 

EOUOOBO  ) 

m 

RAOAk  has 

CHANGED 

TRACK 

c 

EOOOOOA  ) 

m 

( flag  not  USED 

) 

V 

c 

c 

E0200BO  ) 

m 

errors  2 

AND  S 

ARE  BOTH  present 

V. 

T»BUFF( iStHI  ) 

•TLIFT 

0T»BUFn  IG,RC)-BUFFUS»MI  ) 

meuFFis 

»R1)#GT.I,) 

ITR*  IM00EI2)  ti  skin 

IFCpUFF  ) tGTtSt  ) I TR«  IMoOEm 

IF  iBUf  F (^,<41  ) tUTtOt  ) I TR«  IHODE  < 3) 
t TQ( I )*BUFF (6  ,<U  ) 

lTD(2)*eOOLCt)UFF(7,4(i  ) ) 

I TD(3  )aBQOL<BUPF  ) ) 

RaBUFFI 1 |R| ) 

AoBUFFI3»Hl )•  RTDO 
E-BuPFIRfRl )•  RTOO 
0DR5I  • 0UFF|9,H|) 

0EL2(  ■ Burn  lOiHl  ) 

SGDR2l«  Burn  I 1 |R1  ) 
SVeL0.UBUPFn2|Rl  ) 

OR2)  aeurnutRi) 

VEL«HOlDU  i2il  ) 

0EV2I  ■HOLOn»l»3)  • hold:! 

DEVS)  •HOL0(2ilt3) 

INPTS»|NPTS*I 

|F(BU^riSiHntLT*Ot  ) GO  TO  3010 


«2»1  } 
HOLO(1»2*|) 


ti  BEACON  TrK 
NO  TRK 
Q VEL  HOoE 


A)  range  fTi 
A»  A2  o^G 
>9  EL  oCG 
9 range  ACC 
® VEL  ACC 
« RVEL  iSl6 
9 OOP  IsU 
9 range  vCL 

9 ROOT-  RAP  VEL 
9 ROOT-  rAB  VtL 


9 ONLY  COUNTERRqRS  RaDAR  INTRACk 


test  COUNT  FOR  GOOO  OOPP  EITHER  APPROACHING  OR  LEAVING  RADAR 
COUNT  HA'.’E  only  FOR  VEL  G*Tt  ABSUOO  FP5) 


S5T 


|F(ARScH0LU<2i1 »3) ) tLTt 100*0)  GO  TO  S&6 
inHOL0(2ili3)*LTi-|OO»O)  GO  TO  hb7 
IPOSPTalPOSPT*) 

|F(  A6SI0EVS1  ) *LT*  l<i  *3l)  I VPOSaj  VP0S9] 

GO  TO  S* 

CONTINUE 

INCGPTalNEGPT*) 

IF i ABStOCVSl ) tLT  » )H t 3) ) I VNEG* f VNE G«  1 
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9 IF  transfer  circling 

9 Target  receeding 

V TAt.GEr  RECEEQING 

9 Target  receeding 
9 Target  approaching 
tf  Target  approaching 
9 Target  approaching 


S56  CONTINUE 
C 

c 

IF ( ITD(2)  tEQ* I > ISPTS-lSPTS+1  « COUNT  DoF  Sk I N FTS 

IK ( I TD< 3)  tEUi 1 i IBPTSbIBPTS*!  Q COUNT  OQP  CQMO  PT3 

IF  (BUFF  ( 1 I iHU  •LT«6,00)  IRlilO  "IRblG*! 

IF  (BUFF  ( 12  »**M  «LT.  U2U)  I OS  1 G >10S  IG*I 

C TtST  ERKOR  MHEN  DATA  IS  BAO  BUT  OvCS  VALID  SAYS  ITS  6000 
C 

lF(ABS(DCV&l).GEtlH,3l.ANDtlTD(|l«CQil)  lOBFIalOBFI^I 
C 

C TEST  ERROR  MhEN  DaTA  IS  GOOD  BUT  oVES  VALID  SAYS  ITS  HAD 
C 

lF(AB5(DEVSll«LT.lH.3ltAND*lTDfl<'*CQiO)  1DGFO«IUGFO*I 
3010  continue 
C 

IF ( IOPT»NE  • 1 > GO  To  399  « SKIP  PRINT  ERROR  ONLY  LOGIC 

DO  3S0  JJ«1  i2h 
J«I *JJ 

DO  3SI  I«1  I IH 
3S1  5PR I ( I , JJ  >«SPRl ( I iJ t 
3bO  CONTINUE 

5PR  I { 1 ,2S) «BOOL  { lERR) 

SPRI (2i2S>*T 

SPRI  (3i251*iOT 

SPRI  Mt2S)«B00U(  I R| 

SPRI  IS,2bl>OOOLl  WD(  1 > I 
SPRI  (6,2SI"B00L( ITD|2)  I 
SPRI  ( 7 ,25)*POOLn  TOO  I ) 

SPRI  ( 6,2&>*P 

SPR  I (9,2S) MA 
SPRl  Ud  i2S  »-E 
SPRI  ( 1 1 i25)«00»51 
SPRI  I I2i2S  >*0EL2 I 
SPRI  ( I 3i2S)«SGDR2I 
SPRI  ( m OS  1 ■SVEL2  I 
SPRI  USf2SM0R21 
SPRl(Ui2S)>VEL 
5PRl(I7i2SlaDEv2I 
SPRI  ( le  i2S  XDEvSl 


IKP»IKP*1 
1 TEST«0 
DO  3S2  1>1,2S 

ITeST«iTEST  * (BOOL (SPRI C I , I n ) 
]F(  ITESTiEQ>0)  go  To  398 
|KP«  0 

ITEST«lTEST»(BOnL<SPRI ( 1 i2b)  I) 
IF(  IVEST.NEOI  go  To  399 
DO  3S3  !•!  i2S 
ICRRB  bOOL  (SPRI  Mill 
T ■ (SPRH2»i 

OT  ■ (SPRI (3,1 

ITR  s BOOlISPRI ( R I I 
|TOU)*BOOL'SPWMSiI 
ITD(;>)aB00L<SPRl  (B,* 

1TD(  laBOOLlSPRI  ( T,  { I 


9 NO  errors 


I9IF  TRANSFEH,PREVt  2S  AL  RE  AO  Y DuMPE  o 


<9  CUhP 
Q 

<i  PREVIOUS 
9 

9 2S 
Q 
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R m 

5PNI ISiI ) 

.2 

DATA 

A B 

SPNI 1 9» I I 

9 

E B 

SPRM  10,1) 

9 

Points 

OORSJ  a 

SPRI M I , I » 

0Eu2l  • 

SPRI 112,1) 

5G0R21* 

SPRI 1 13, t ) 

SVEL21* 

SP«H  IN,  11 

0R2I  • 

SPRI (15,1) 

v£l  • 

SPNM  )*,/  1 

DEv2l  ■ 

SPHH  17,1) 

DEvbl  B 

SPRM  13,  I > 

IF(  It'RK  ( 

>EQ«0)  GO  TO  3S5 

ri'U  t20CO>  ICRft  >T»DT  » I TV  , f rO  f W t r ^0  ( 2 > i 1 TO  t 3 > IN  » Aie  lOOKbl  iDE<-2| 
U$G0N2l«Svet.2((0i«2UVtLi0Cv&lt0CV2J 
GO  TO  3&A 
3&b  CONT  tNuE 

fi'^\ft^bt2D0\  ) (/MH,  r lOTi  t TH  I I TP  I 1 > I f TOf  2>  flTOM  j t 1 A 1 1 1 DORS  1 ioEU2l 

i f SGPR2{  |SVL(.2l  tON2l|Ve.LiOCVSl,PeV21 
ibb  CONTINUE 

3S3  continue 
SO  TO  39fl 

continue 

|PI  IERN«EQ«0)  go  to  naq 

WRITE  lGf2DOC)  KNR  ,Ti0riITKiIT0Cl)ilT0(2)ilTD(3)  iKi A lEf DDR&l iOEL2l 
1 rSGDN2i  |SVEU2]  i DR2  j « |P€  VS  I |0EV21 
GO  TO  NAl 
N60  CON1 tNul 

WRITE  t’A  |2G0  1 f [NaM  iTipTt  1TN|  I TP  ( 1 > f I TP  < 2 > i I TP  f 3 / iRi  A lE  iDPRSl  iDCl2l 
1 f 5GPN21 i5V£L2l  iOR2l  iVCL  |DEVSt  , PE V21 
R*I  CONTINUE 

if  8 continue 

IF ( IOPT*ME* II  GO  To  SOI 
IFC IKPiUTi 100)  GO  TO  bOO 
WRITE(6t20Cl I 

WRITE  (6f20pn  inaHi  r ipr  1 1 thi  I rp(  1 1 1 jTpui  1 1 rooi  iRiA  lEiPORS)  (OCl;) 

I f 5GOR2I f SVEE21 . 0r2 1 , VEt , DC VS { , OE v2 l 
aR1TCI6|200J ) 

501  IKP-O 

5c.-;  continue 

c 

IFIBUFMStNl  I iLTiD#  )Q0  TO  300 

IPTCwiPrC*)  H TOT  NO  RTS  COUNTER 

IF  ( IcRR.NE.O ) IRTCE  • IPTCE  ♦]  9 TOT  CRRORS  COUNTER 

IF  ( It f 31 •NL»0) JRTC0*  IPTCe  ♦!  9 BIAS  ERROR  cOUNTiff 

300  continue 

IFfCNOFILtEGtl ) GO  TO  R62 
C 

IFUTIHE»EU*Oi  ANOiTSToRfET*  1 <01  GO  TQ  tOO  9 1 F TRaNSFR  CQNT«RR0C* 
tF(inME*EU*0)  GO  TO  nh3 

I T5T«iriMC»20  ra  process  all  RCO.  ITIme  points 

I00P2«|00P2*1 

|FU00P2*LT*  ITBT  ) GO  TO  lOU 
GO  TO  RR9 

RH3  iFIfiUFF nSiRI J .lT.  TSTOPI  GO  TO  100  Q PROCESS  TO  TSTOP 
60  TO  N*(9 


continue 
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MH9  TNUM  ■ lPTCr*lDO 
FDEN  • IPTC 
ERPH  ■ FNUM/FDEN 
FnUM  ■ 1PTCB*»00 
VLPR  « FNUN/FOEN 
«RlTE(!i,2n) 

WHlTEU|21b>  INAM 
WRITCt6,212> 

ir(tNNN«Ut*&nARITE  <6|H&OOI 
H^OQ  FOPHAT(»SH  RADAR  NEVtR  SHOrCO 
TF  « BUFF  US  »R  I ) 
ttRneUtRSOl)  TltTF 
HSOl  F()RHAT(lRK  START  OF  RUN"tFlOi 


Q calculate 
ra  percentages  oF 
<9  total  errors  and 

r;  BUS  ERRORS 
rf  ANO  PRINT 


beacon  or  SKJN  track  - NO  DATA  ) 


lUH  |End  or  run«ifioo) 


ARUt  U»2U) 
wR1TEU»h602)  IPTC 

HS02  FORMAT137H  NOi  OF  PqINTS  OF  UiTRACK  RAQAR  OATA«iIiO) 

HS03  F0RMAT^37m  TOTAL  NUMBER  OF  RAO*  PTSi  PROCESSED*  i U 0 » 

WRITE(G>3007I 
WKlTEIGiRSOS) INPTS 
WHITCUi2I2) 

ARlTEUiHSO)  ERPRiVtPR 
WR1TE(6i212) 

Rbo  formati22h  Percentage  of  ehrors«,fs, i iJnh  percentage  of  bia 

ISEO  DAtA**FSU> 

WR1TE(4i3007> 

FNUM«  lOBFUlCO 
EOBF  1*FNUM<'FDEN 
FNUM«  lOGFO*lOO 
edgfo»fnum/fden 

soil  F0RMAT(H2H  percentage  of  invalid  OVES  a/good  data  ■ *FS*I> 

3012  F0RMATIN2H  PERCENTAGE  OF  VALID  DV£S  A/BAD  DATA  ■ iFStll 
wRiTElb, 30111  EOGFO 
AR]TEUi3007) 

WRITEUi3012>  ED9F) 

WR1TEU|2I2) 

ARlTEUt3007l 

SDREL  b(lDSlG*lOOQ)  /IPTC 

SDREL  •(5DREL/10*0)**OB 

5RREL  bURSlG«1000)  /IPTC 

SRREl  blSRREL/lO*OU*OB 

ARlTE(A|>4B0e)SRRCL 

ARITE(6i3007) 

ARITE(4|HB09)SDREL 

HB08  F3RMATU2H  RELIABILITY  OF  LOW  NOISE  ON  RANGE  SERVO-  .FSiliBH  PERCe 
INT  > 

HSQ9  F0RMAT(H2H  reliability  OF  LOA  NOISE  ON  DQPP.  SCRVq-  iFSUiBH  PERCl 
INT  ) 

AR|TE(4i212) 

GDFLG-  I00*n**(EDQF0*C0BFl  I 

301H  F0RHAT(31H  RELIABILITY  OF  VALID  TRK  FLG«  »rS«I|9H  PERCENT.) 
ARlTEUi30lN)GDFLG 
AR|Te(«i3007) 

XPER*USPT5*10001/ IPTC 

xpcr«xper/iq*o  *>US 

AR1TE(6iHS0S)XPER 

HSOB  F0RHAT139H  RELIABILITY  OF  OOPP  SKIN  RETURN  FLaG-  ,F5*|»SM  PERCENT) 
AR  I TE (6i3007 ) 
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XP£R"< IBPTS*JOOO)/|PTC 
XPER"XPER/10«0  4>«0S 
WRITE(6»M5Q6)  XPER 

*♦506  F0RMATO9H  RELIABILITY  OF  DOPP  COhO  RETURN  FLAfi*  PERCENT » 

RRlTEUt212) 

IF( IPOSPT«EQ«0)  60  TO  5B6 
VPP«< IVPOS»IOOOI/IPOSPT 
VPP»lVPP/10t0U*05 
RR  I Te (6iB70) VPP 
NRITE(6i3007) 

570  F0RHaT(53H  reliability  OF  UOPP  VEL  AHEN  TARGET  IS  RECEEOlNG  ■ iFg 
1*1 »8h  PERCENT) 

558  IFI  INEGPT»E0«0)  GO  TO  5S9 
VPP*I IvNEG*IOOO)/1NE6PT 
VPP«(VPP/I0iO)*f05 
WR1TE(6*&71 )VPP 

571  F0RHaT(63H  reliability  OF  DOPP  VEt  WHEN  TARGET  IS  APPROACHING*  iF& 
ItliBH  PERCENT) 

559  continue 
MrRlTe(6t2I2) 

XPER«I  IPTC*I000»/INPTS 
XPER^XPCR/lO**)  *»05 
NRITEUiRSOR)  XPER 

MSOH  F0RHAT(29H  radar  TRACKING  RELIABILITY*  ,FBfli6H  PERCENT) 
WRITEUi3007) 

GOATA*  100*0«VLPR 

3013  F0RHAT(3IH  reliability  of  VELOCITY  DATA*  iFS*1i9h  PERCENT*) 
PR1TEI6|3013)G0ATA 
RRITE(A»213) 

RR|TEU»212) 

WRirC(A|212) 

212  FORMAT!//) 

213  FORMATIlHl) 

219  F0RHATI27H  input  DATA  CARD  FOR  CHECKS) 

215  F0RHAT(29H  radar  TO  BE  SEARCHED  • »A9) 

21*  F0RMATI29H  TIME  TO  START  TAPE  ■ iFlD*3> 

217  F0RMAT(2HN  time  TO  STOP  RUN  • |F10«3) 

218  F0RMATI29H  OUTPUT  PRINT  INTERVAL*  iF10*3) 

219  F0RMAT(29H  PRINT  TIMES  BIASED  BY*  tFlQo) 

220  F0RMAT(**H  when  THE  FOLLOW  OPTION  IS  USEDi  iT  WILL  OVERIDC  TSfART 
1 AND  TSTOP* ) 

221  F0RHAT116H  PRINT  ONLY  THE  iI2|22H  < 1 ST i 2ND i 3R0 i 9 Th * E Tc ) i 
1I3H  interval  of  tllOilBN  SECONDS  OF  DATA*) 

222  FORMATlSOH  ONLY  FLAGOED  ERRORS  WILL  BE  PRInTEo  WITH  lOPT  *1  ) 

223  FORMaT(*SH  all  data  WITH  AND  WITHOUT  ERRORS  ARE  To  BE  PRINTED  WITh 

1 lOPT  * 0) 

229  rORH.'.T(39H  ERROR  IN  lOPT  COL*  BOiSHQULD  BE  1 OR  0) 

225  FORMAT!  I 30H  t ••••••• 1 •••••••••• 1 | • 1 1 1 • 1 • I •••••••• t ••••••  I *• 

I I 

2»******»I  ) 

3050  F0RMATM30H  I ••  t ••••  I ••••••••••  1 •••••  »t  | R 0 5 B I *•••••••  I •••*••  I •• 

I 

I I 

22*  FORHATdSOH  lERR  FlgI  TIME  J O^L  T I TRACK  1 RANGE  1 AZ  1 

lEL  I DDRNG  1 DVEL  1 SGDR2 1 1 SGVL2 1 I RNGDOTSll  RAWyEL  III  VELeRR  1 

2 VELERR  I ) 

227  F0RMATII30H  lOR  NAM  I !SEC)  I (SEC)  I FLAGS  1 |FT)  1 !0EG)1  ! 
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IDFGH  K/S2  I r/S2  I TPS  I FP5  I F.PtSt  I F»P.S.  Ill  5lPTS»  1 
2 2 IPTS.  I J 

2000  f0R^AT(2M  t|l6i2H)  ,F7«i«lH  ,Ali3(th  ilpilh  iFS^Ot 

I7UH  ir6t2)i2<lH  iF7<n)t2llH  ,F6*2),2Uh  iF°*2>t2H  i2(lH  tFS,2)) 
2Q01  FOHKATtSH  »F70»JH  iAI,3Uh  ,FB*Oi 

17UH  iF6»2)i2tlH  iF7»0ii2WH  ,F«*2)i2UM  iF9f2)i2H  i2UH  •Fei2>) 
2002  FOKHAT(AH,|6,<4F10t3,t|OtUU(19»m 
STOP 
END 


COMPILATION!  2 DIAONOSTICS. 
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1 


I 


LOGRB 

•0l/<3/7N*i  M591  j7  C,o) 


] 

r 

j 

r 


LLL 

ENTRY  POINT  00lD*S 

Dt  COOEt 

> ) 0011S3I 

DaTAIQI  00332HI 

BLANX  COMMON  121  OOQOOO^ 

f ERENCES 

(BLOCK,  NAME) 

R38 

1 QNMCNT 

(BLocxi  type,  relative  lqc at IOJP^^ME  1 

OHO 

tOOL 

bOOl 

000032 

1306  . 

0001 

QOOOS* 

1376 

0001  ' 

'OQOIO* 

1636 

113 

1*2* 

DODl 

OOOIH* 

17H6/^ 

0001 

OOOIBO 

(776 

0001 

000)63 

2036 

210 

2206 

0001 

0002(3 

22^ 

ocoi 

0002H* 

23*6 

0001  ■ 

000260 

29(6 

320 

2*7(, 

0001 

000330/ 

'17T6 

0001 

00027* 

300L 

0001 

0003*3 

3026 

S7H 

30Hl 

0001 

OOOXl 

30*L 

0001 

000997 

3296 

0001 

‘ 000967 

3276 

S7* 

3*U 

QOOl 

OB<*0* 

3*96 

0001 

000*91 

3796 

0001 

00072* 

91*6 

7*1 

H306 

0001 

^/000025 

SOOl 

0001 

000197 

601l‘ 

'Wor 

eroo202 

6Q2l 

03S 

00^ 

D OQOn* 

A 

0000  D 

OOQl 10 

CON 

OOQO 

0 0010*9 

0 

SIS 

DVV 

0 00300* 

ovvv 

0000  0 

OQOOOO' 

02l 

oooo 

0 000022 

031 

0** 

OBl 

^^0000 

1 C03220 

1 

0000  1 

003216 

|PLA6 

0000 

0032*2 

INJP: 

221 

J 

0000 

1 C03222 

K 

0000  1 

003M7 

kset 

croQiO 

R 003223 

SUM 

207 

t\yt 

0000 

R 00320H 

XPT 

0000  0 

003211 

XTlHE 

0000 

R 003229 

ZSTX 

i: 

! 

i: 

i: 

c 

E 


subroutine  lLL(NPO|NT,DeRKET,rtNVARfXRAll,XSi'.SIGXtD>(SHtDSlGX, 

1 UOXSH  |DOS IGX iCVXR ) 

dimension  D2ll3i)}tD21(3i3)tOslt3iSt(DBIfli3) 

dimension  CONI3)|A{Bl,3l,OI3,8ntDV(.>,Sn,OVv<DOI)tDVVV(3|2p 

dimension  r INVRRtBl ) iXPT (3t 

DOuBUE  PRECISION  CON t X 1 NC i XT  I HE i A , 0 i D V 1 0 VV t DVV V i T 1 ME2 

DOUBLE  PRECISION  02  I 1 1>3  1 iDB  1 t DB  1 

DATA  d2I'^0*7i32NSSaMS0N|H6MD*01  > -0 • I IBBSBmS? aq  1 0]  AID* ) * « 

A -0>BS380BMa>4l28Bi)0S0*00i*0*lBBS3BMB7A01QlA8t}-I6i 
B O»Sl9N8OSl9H8OS19BD*0OtOi3l0>HSBl2|03*72HO-l7i 
C •0t6S3BDBH3Nt2B8oo60«00tQf3|Q)HB8I2l0S*72M0*l7'i 
D 0«t07S6|NB7«tlMlB7AD*00/ 

data  DSWOi  lOI  122hS929B2102d*01  |0*382  I5439S&M7Q7hB0-|7  , 

A 'Q>2022NH9i8S9A'«2030*OOiOi3821&4  39SSh7o7mBD*17  I 
B Ot  1*12903 2 2S80ANS30*00(-0*  1 30  1 9622m  199  1 733.0*  I 7 i 
Co0i2022«tN918&9AM203D*Q0t-0*  13019S22h199  17330-17. 

0 0»727070N8235HI30Bd-01/ 

OAT*  061/0*B3*32872B97S1&OSD-01 i0**7222Si%3B«2379H0-|8i 
A -0*HS30I 139323*S3920-q1 iO t *7222tM&3S*2379MD- 1 B , 

B 0*3«l9909S022*2*tMBD*0l  i-0*3*H  1221  20*79*2:60- iBt 
C *0i9S301  139  323*539  20-01  i-0«3*9l2212Q*79*2t'&D'-lBi 
D 0*HM1HE9*Q2709Q!373D-01/ 

OAT*  OBl/0*82**0B79BHl&S927D-O2fO«9222M3*7*970QSolD-l8« 

A -Ot 1 129»  >HSB70*B3B*D-o1 |0»922293*7*9700I010-1B| 
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1 


\ 


r 


0 0«9033M  236675700260-02,-0  *67  6 I 1 93H‘*?7257H  20- I 8 , 

C -Otl  129  69H56  7D6B36HD-DI  * -0  • 676 1 1 93i<«l97  257  Hi  0- I 8 , 

Q 0»2778H8383689H6900-0|/ 

DATA  irUAQ/  1 / 

c 

c 

C SYmBOuIC  DICTIONAKY 

c 

C NPOINT  • NO*  OF  PTS*  OVC»  PHJCm  CURVE  IS  TO  BE  FITTEO* 

c derkCy  • 1*0,  cohpute;  derivatives  oF  second  oeq*  eq> 

c derkey  ■ 0*0,  NO  derivatives  computed,  sjhply  fit  Input  data. 
c FINVAR  ■ input  data  ARRAY  WHICH  IS  ALWAYS  8I-PTS. 

C XRAW  ■ RAW  INPUT  VALUE  *T  THE  Hl-PT* 

C XSM  ■ SMOOTH  VALUE  AT  THE  HI»PT. 

C SIGX  ■ variance  of  position  at  The  Rl-PT. 

C DXSM  • smooth  FIRST  DERIVATIVE  AT  T^^E  Hl-pT. 

C DSIftX  • variance  of  first  derivative  at  The  MI-PT. 

C ODXSM  • smooth  second  derivative  AT  THE  **l-PT. 

C ODSlGx  « variance  QE  SECOND  DERIVATIVE  AT  THE  hUPT* 

c CVAR  • variance  of  SECOND  deq.  curvE  fit. 

C XINC  • TIME  incrementing  CONSTANT  FOR-- A--M A T R I X , 

C XTIME  • TIME  generated  FOR  — A--M  aT  R I X . 

C TIME2  • TIME  SQUARED  GENERATED  FOR  — A--M  aTR  1 X , 

c Atiij)  ■ transformation  matrix  to  transform  The  constants  of 

C A SECOND  degree  EQ.  TO  SMOOTH  POSITION. 

C C0N(3)  > Atl,3)  • 1.0,  or  transformation  matrix* 

C C0NI2)  ■ AM, 2)  • time  of  PT.  RELATIVE  TO  STARTING 

C time  of  curve  fit. 

c coN(i)  > Ati,i)  • time  squared  or  pt.  relative  to  starting 
C TIME  OF  curve  fit. 

C lx  ■ NO.  OF  PTS.  second  DEGREE  EU.  IS  TO  BC  F|TTeD  OvER. 

C KSET  ■ NO.  OF  PTS.  TRAVERSED  TO  THE  RIGHT  AND  lEfT  ABOUT 

C MIO-PT.  REFERENCE  TIme*  REFERENCE  TIME  STARTS  AT  ZERO. 

c This  symbol  «-)  is  not  a minus  sign»  it  is  a dash  in  the  eqs. 

c D21  ■ (A-TRANSP0SE*A»-1NVERSE  matrix  FOR  A 2l-PT,  CURVE  FIT# 

C 031  B (A-TRANSP0SC«A)-1nVERSE  matrix  for  a 31-pT,  CURVE  FIT# 

C D51  » U-TRANSP0SE*A).1nVERSE  matrix  FOR  A SI-PT,  CURVE  FIT# 

C 081  ■ U-TRANSP0SE*A)-INVER5E  matrix  for  a BI-pT,  CURVE  FIT# 

c d ■ Matrix  to  transform  si-pts.  of  data  into  the  3 parameters 

C OF  the  second  deg.  EQ# 

C DV  « matrix  To  transform  Sl-pTS.  OF  DATA  INTO  THe  3 PARAMETERS 

C OF  The  second  DEG#  EG. 

c ovv  > matrix  to  transform  31-pts.  Of  data  into  ThE  3 parameters 

C OF  THE  second  deg.  EQ. 

C DWV  • matrix  to  transform  21-PTS.  of  data  Into  the  3 parameters 

c of  the  second  deg.  EQ. 

c 
c 
c 

c The  set  of  INSTRUCTIONS  BETWEEN  STATEMENT  NUMBERS  500  AND  300 

C ARE  EXECUTED  ONLY  ONCE.  INITIALLY,  IFLAG  iS  SE^  To  ONE  IN 

C ORDER  TO  EXECUTE  SUCH  INSTRUCTIONS. 

C IX--1S  SET  FOR  an  81-PT.  CURyC  FIT. 

c ;cset--is  set  to  the  lower  and  upper  bound  span  abOf\jt  mio-pt. 

c reference,  reference  time  starts  at  zero  and  increments  by 

C 0.05  SECONDS. 

C SET  C0N(3)  TO  ONE  FOR-- A — M ATR  1 X . 
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c 


lF-(  lFuAG«E.QtO)GO  TO  3O0 
I X > 81 
KStT  ■ Ml 
CON  I 3 1 ■ I .0 

c 

C generate  time  iNCKEMEnTlXlNC  » I TlMEUONlJl).  ANU  TIME  SQUARED 

C (CON  (II)  FOR  — A ( 1 , J)»-HATRU  • 

C LOAD  computed  VAUUESICCN( 1 I)  I NT 0*- A ( 1 , J ) .-M A TR I X i 

c 

500  DO  5 I>)  t I X 

XINC  ■ I-KSET 
XTIME  • 5tOD-u2«XlNC 
Cun ( 1 I ■ XT  I ME • *2 
CON  I 2 I ■ XT IME 
DO  5 JBI »3 

5 A( 1 , J)  • Coni j) 

C 

c Select  the  si-pt,,  or  3i-pt*i  or  2i-pt.i  path  determined  bt 

c The  preset  value  of  •-ixi 

c 

IF( tXiEQ«5l ) GO  TO  SQ 1 

1F( IX«EQ«31  } GO  TO  S02 

ir( IXtLQt2l ) GO  TO  503 

C 

C generate  the  DATA  TRANSFORMATION  HATR 1 X^-Q--FOK  THE  61-PT*  FiT* 

C 

DO  7 1>1  ,3 
DO  7 jaiiBl 
D( t , J)  ■ 0.0 
DO  7 X« I i3 

7 D(I,J)  > D ( 1 , J)  *Dm  I ( I |K ) *A( JiK ) 

IX  « 51 
KSET  « 26 
GO  TO  500 
C 

C GENERATE  THE  DATA  TRANSFORMATION  M A TR I X --0 V- -FO H ThE  51-PT*  FIT. 

C 

50:  DO  fl  !■! .3 

Do  0 JB  1 (51 
DV(  I lU)  B 0(0 
DO  a KB  1 (3 

a Dv(!iU)  • DV ( 1 I J > *D5 1 ( 1 (K I • A ( J |K ) 

U • 31 
KSET  B 16 
GO  TO  bOO 
C 

C GENERATE  the  DATA  TRANSFORMATION  M A T R I X- 0 V V “FOR  31-PT.  CURvE  FIT, 

C 

502  DO  7 ia| ,3 
DO  9 JB  I ,31 
DVV ( I , J ) B 0,u 
DO  9 KB  1 (3 

9 DVV(1,J)  B DvV n I J t «D31 < 1 iK ) tA ( J,K ) 

IX  B 21 
KSET  B II 
GO  TO  500 
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c 

c 


generate  data  transformation  MATR  I X-OVVV-rOR  21-f-T,  CURVE  riT 


Bo;)  DO  10  !■!  »3 
DO  10  t2  I 
Ovvv({(w)  • 0*0 
00  10  K«  I .3 

10  DVVV(liJ)  ■ 0VVVU|J»*021(l|K>*A<J|l<» 

c 

C SET  IFLaG  to  zero  TO  SKIP  STATEMENT  NUMBERS  500  THROUGH  3o0« 

C 

IFL*G«0 

c 

C SET--SUM--TO  ZERO,  (SUMMATION  OF  RESIDUALS  SQUARED) 

C SET  PARAMtTERS--AiBiC-(XPTll) 1--0F  SECOND  DEGt  EQt  Tq  2ERQ» 

C 

30C  SUM  ■ 0,0 

XPT ( 1 ) " 0,0 
XPT(2>  ■ 0.0 
XPT(3)  ■ 0.0 
C 

c Select  the  appropriate  path  for  the  no.  of  pts,  usep"  in  curve  Fit 

c 

JF(NP0INT*EQ.R1>00  to  306 
.«‘(NP0InT,EQ,S1  )G0  TO  30R 
IF  (nPoINT,E«, 31 ) 60  TO  302 

c 

c COMPUTE  Parameters*  a,biC-(xpt( i i ) for  second  deg*  eq. 

C FOR  A 21-PT.  CURVE  FIT, 

C LOAD  SMOOTH  value  AV  Rl-^PT,  INT0*-XSK. 

C COMPUTE  SUMMATION  OF  RES1UU*'S  SQU AR EO--S UM • 

C 

DO  301  1*1 i3 
DO  30)  J*31 ,51 
K«J-30 

301  XpTtl)  ■ XPT ( 1 ) ♦OVVV ( 1 ,K ) *F invar (J) 

XSH  » XPT-3) 

00  hOi  !■! ,21 
X NC  • 1-1 1 
XTIME  ■ d,0D-02»XINC 
T1ME2  ■ XT1ME#»2 

2STX  a XPT ( I ) «T 1ME2*XPT C2 )*XTIME*XPT ( 3) 

K«l+30 

*401  SUM  • SUM*(FINVAR(K)-2STX)*«2 

c 

«.  LOAD  CURVE  VARIANCE  INTO  — CVAR. 

C calculate  "ID-PT,  variance  of  POSITION— SlQX, 

c 

CVAR  ■ 5UM/1S,0 
SIGX  1 0»597508263E-n2«£JM 
C 

C lF-*Ti,  RKEY--EOUALS  ONE,  LOAD  IN  FIRS)  DER  1 V AT  | v E»  ( OXSM  1 , 

C C0NPU1E  variances  QF  FIRST  D£R 1 V AT  1 V£ - ( DS I GX  ) , AND  VARIANCE  OF 

c The  Second  orRivATivE-iODsiax), 

c 

STIC*  TH*i  TEST  FOR  EQUALITY  BETREEN  NON-INTEGERS  MAY  NOT  Be  MEANINGFUL, 
ir(DCRKEY,EQ.Qf OlGO  To  100 
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OXSM  • XPT(2) 

USUx  ■ 0«2886C>C288E-01»SUM 
DOSIGx  ■ 0. l&BM99oj3E*0J*bUM 
60  TO  99 


C 

C 

C 

C 

C 

C 

302 

303 


H02 

C 

c 

c 

c 


c 

c 

c 

c 

STIC* 


compute  parahEtehSi  a »b iC- I xpt I n > For  second  oCQt  Eq* 
r OH  A 3 1 ^PT  ♦ CURVE  FIT* 

UOaO  smooth  value  at  *tl“PT»  1NT0-*XSM« 
compute  summation  of  residuals  syuAREO— Sum. 

DO  303  l«Ii3 
DO  303  JM26|S6 
K«J-2b 

XPTll)  « XPTt  I I+DWM  |K»*F  INvARljI 
XSM  ■ XPT(3> 

DO  902  tPl  t31 

AlNC  • 1^16 

XTIMg  • S.0D-02*XINC 

TIME2  * XT1HE**2 

2STX  m xPT ( 1 ) *T  IME2*XPT t 2 » * XT  I ME*APT t 3 ) 

K*I *25 

SUM  • SUM* I F I N V ah ( K I *Z5T X I • »2 
LOAD  CURVE  VARIANCE  INTO  — CVaR* 

calculate  MID-PT#  variance  Of  P051T10N--5IGX* 

CVAR  » SUH/20.0 

SIGX  • 0*259fr6802SE-02»SUM 

IF  — OERKEV— E'JUALS  one,  load  IN  FUST  OER  1 V AT  I V £- ( 0 »SM  ) , 
COMPUTE  VARIANCES  0^  FUST  D£R  I V A T I V E - ( 0 5 I Q X I i AnD  VARIANCE  OF 
THE  SECOND  OE R 1 V A T 1 V E - 1 DOS  1 <»X  ) • 

THE  test  for  EQUALITY  9ET*lEEN  NOn-InTEGERS  may  NOT  BE  MEaMNQFuL. 
IFIDERKEY.EG.O.OUo  to  100 
DXSM  a XPTt2l 

DSIGX  • 0.57603686nE-02*SuM 

DDSIGX  • 0. I H996069flE*0J*5UH 
GO  TO  99 


C COMPUTE  PARAMETERS,  A,Bi  -'XPTIll)  FOR  SECOND  dEg*  EQ. 

C FOR  A Sl-PT.  CURVE  t T # 

C LOAD  SMOOTH  VALUE  9J-PT.  lNTO--XSMt 

C COMPUTE  summation  RESIDUAlS  SQU aRCO“»SuM • 

c 

309  00  30b  l«Ii3 

DO  30b  vlPI6,66 


K«J-Ib 

305  XPTU)  • XpT  I n*OvM  ,K  >«F  UvaRI  J» 

XSM  • XPT(3» 

DO  90s  UliSI 

XINC  ■ U26 

XTIHE  • S«0D-02**INC 

TIMC2  • XTiMt*»2 

ZSTX  a xPT ( I 1 *T  |ME2*XPT (2 »*XT IME^APY U) 

X*  l*lb 

<i03  SUM  a SuM*UUV*RtK»-ZSTxU*2 
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c 

C tOAD  curve  variance  INTO--CVARt 

c Calculate  mip-pt«  variance  Of  positiOn--sigx« 

c 

CVAH  « SUH/HBtO 
SlGX  m 0t9  I970>**98E»03*SUM 
C 

C IF— .OtRKEY^-EflUALS  ONE  i UOAO  IN  FIRST  DER  I V AT  | VE- I ()XSM  I , 

c Compute  variances  of  first  derivative-iosigx) , and  variance  of 

C The  second  DERI  vat Ivt-(D0S1GX ) t 

c 

STIC*  the  test  for  equality  between  NOn-INTEGCRS  may  not  Be  meaningful* 
IF(DERKEY.£q*0*0)G0  to  lOQ 
DXSH  ■ XPT  (2) 

DSIGX  • 0*7SR  lR7eoHE«03*SuM 
ODSIGx  ■ 0«69699060<4E-02*5UM 
GO  TO  99 

c 

C COMPUTE  PARAMETE'^Si  AtBiC*UpT(|n  FOR  SECOND  DEg*  EQ* 

C FOR  A ai-PT.  CURVE  FIT* 

C LOAD  SMOOTH  value  AT  *4l-PT*  1NTO--XSM, 

C COMPUTE  summation  OF  RESIDUALS  SQUaRED-*Sum* 

c 

306  DO  307  I«1 »3 
DO  307  J»l  fSl 

307  XPTUJ  • XPTI  U*DU  ,JI*FINVARIJ) 

XSH  • XPTI3) 

DO  HOm  I»1  |BI 

XINC  • 1*41 

XTI.ME  ■ 5*00-0?*XInC 

TI’(E2  ■ XTIMEa*2 

ZSTX  « XPTU  )«T1ME2*XPTU)«XTIME*XPTI3> 

ROM  SUM  • SUH*(F1NVARI 1 )-ZSTX>««2 

C 

C LOAD  CURVE  VARIANCE  INTO««CvaR* 

C calculate  MI0»PT*  variance  of  POSIT10N-«SIOX* 

C 

CVAr  m SUm/76.0 
SIGX  m 0*3S62 IS07 iE«03«SuM 
C 

C IF--DERKET--EQUAL5  ONE*  LOAD  IN  FIRST  OER  I V AT  I VE- ( D»^SM  » , 

C compute  variances  OE  first  DERIVATIVC-IDSIGXI  , AND  VARIANCE  OF 

C THF  second  DEr I VAT  I VE- (DDSIQX ) • 

C 

>•  . . : tuUALlT’.  0-  .J-INTE6ERS  MAY  NOT  B£  MEANINGFUL. 

V %..w  . .EC*O*0IGO  Tv  0 0 
DXSn  • XPT(2) 

DSIGX  ■ 0.1  158  1 3123.  ^03  ‘ ’ 

DDSIGX  ■ O.M239OOMOSE-0.  •SUM 
C 

C calculate  variance  of  second  DERIVATIVE-IDOKSH) . 

C LOAD  IN  Ran  POSITION  AT  n1*Pt*  InTQ**XRA«I. 

C 

99  DOXSM  ■ 2.0«XPTt  I 1 

100  XRAW  > FtNVARlMl  ) 

return 

END 
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POTENTIAL  BENEFITS  OF  COHERENT  C-BAND 
DOPPLER  RANGE-RATE  DATA 

C-Band  tracking  radars  equipped  v,  r;  Coherent  Signal  Processors 
(CSP)  provide  precision  range-rate  dal^  ^ r^rv:  s applications.  The 
following  discussion  will  attempt  to  display  ‘trity  of  range -rate 
data  In  terms  of  Best  Estimate  of  Trajectory  (BE’i  v 'ovement  and 
feature  extraction.  All  of  the  data  shown  were  * .’lected  by  the  AN/ 

FPQ-6  radar  at  Wallops  Island. 

Figure  1 displays  a classic  problem  for  launch  ranges  - trajectory 
determination  during  staging.  The  velocity  vs.  time  trace  In  figure  1 
1$  from  the  BET  solution  of  standard  skin  track  range,  azimuth,  and 
elevation  (RAE)  data  from  a typical  NIKE  launch  at  Wallops  Island, 

The  computer  software  that  was  used  to  process  the  trajectory  Is  a state- 
of-the-art  Kalman  Filter,  similar  in  design  to  many  range  safety  real- 
time filters  used  nationally.  A familiar  technique  of  destroying  the 
memory  of  the  filter  at  known  staging  times  was  utilized  In  this  data 
analysis  to  assist  In  transient  response  at  burn  and  burn  out,  but 
the  velocity  solution  still  exhibits  classic  overshoot  and  subsequent 
undershoot.  The  velocity  trace  In  figure  2 was  made  with  exactly  the 
same  computer  set-up,  except  the  CSP  range-rate  data  was  weighted  in 
the  solution. 

Figures  3,  4,  and  5 are  from  the  L.'acicing  data  of  a super-critical 
designed  fuselage  that  was  dropped  from  an  aircraft.  The  purpose  of  the 
test  was  to  determine  the  drag  coefficient  of  the  body  as  the  velocity 
passed  through  mach  1.  The  test  was  designed  so  as  to  maximize  the 
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velocity  vector  in  the  direction  of  the  radar  range  to  obtain  high 
quality  range-rate  information  from  the  CSP.  Figure  3 and  4 contrast 
the  solution  of  acceleration  vs.  velocity  without  and  with  range-rate 
data.  In  this  case  the  experiment  was  looking  for  a transient  In  the 
trajectory,  and  the  CSP  was  able  to  contribute  greatly  to  the  feature 
extraction.  Figure  5 displays  the  range-rate  residuals  to  the  trajectory 
displayed  In  figure  4.  The  RMS  noise  (skin  track),  with  four  measurements 
edited,  was  8 cm/sec. 

As  was  discussed  in  figures  1 and  2,  the  classic  problem  of  any 
real  time  or  end-point  filter  such  as  the  Kalman,  Is  overshoot,  or 
equivalently,  filter  lag.  This  Is  particularly  evident  when  solving 
for  high  order  terms,  such  as  velocity  and  acceleration,  from  measure- 
ments of  position,  such  as  RAE.  The  filter  must  process  a series  of 
position  measurements  before  the  higher  order  terms  are  observable, 
and  the  result  Is  a time  lag.  One  such  observable  parameter  Is  aircraft 
bank  angle.  The  computer  software  used  for  most  of  this  paper  has  an 
optional  ((ynamic  model  for  aircraft  tracking,  with  bank  angle  and  longitudinal 
acceleration  being  the  highest  order  terms  in  the  model.  Figure  6 dls- 
plfiys  the  X-Y  ground  trace  of  an  aircraft  flight  tracked  by  the  Wallops 
AN/FPQ-6  radar.  The  aircraft  was  equipped  with  a coherent  transponder, 
and  CSP  range-rate  data  was  collected  along  with  the  normal  beacon  track 
RAE  (gross  spectrum)  data.  Figures  7 and  8 contrast  the  BET  bank  angle 
without  and  witn  range-rate.  Notice  the  time  lag  of  approximately  3 
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seconds  In  the  solution  without  range  rate.  Figure  9 displays  the  range- 
rate  residuals,  here  with  an  RMS  of  2.9  cm/sec.  Some  shadowing  of  the 
transponder  may  have  occured  In  this  data,  since  other  segments  of  the 
track  exhibit  2.0  cm/sec  RMS. 

The  purpose  of  the  above  track  was  to  provide  an  altitude  standard 
with  which  to  compare  the  measurements  from  the  NRL  nanosecond  radar 
altimeter  that  was  on  the  aircraft.  Figure  10  displays  the  solution 
for  altitude  above  the  spherlod  as  determined  by  the  altimeter  (stars) 
and  as  determined  by  the  radar  (solid  line).  The  agreement  between  the 
two  systems  is  generally  10-20  cm  over  this  one  minute  span  of  data, 
showing  the  potential  for  both  CSP  tracking  data  and  microwave  altimetry. 

To  round  out  the  spectrum  of  CSP  applications,  a short  span  (20 
seconds)  of  skin  track  range- rate  data  from  a GEOS-II  track  was  Integrated 
to  generate  precise  range  data.  Figure  11  displays  the  orbital  accuracy 
by  comparing  range  measurements  from  a collocated  laser  ranging  system. 

The  X's  represent  the  laser  residuals  to  an  orbit  determined  by  20  seconds 
of  normal  range  data.  The  dotted  circles  represent  the  laser  residuals 
to  an  orbit  determined  by  Integrated  ranges  from  the  same  20  second  time 
span.  The  Integrated  ranges  were  then  smoothed,  and  the  third  trace 
represents  the  laser  residuals  to  the  smoothed  CSP  ranges. 

In  conclusion,  the  above  data  demonstrate  the  utility,  accuracy, 
and  precision  of  CSP  range-rate  measurements.  Furthermore,  the  advantages 
of  precision  range-rate  for  solution  of  high  order  terms,  such  as  velocity  or 
acceleretfee*  has  been  shown. 
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: ABSTRACT 

j 

C-BAND  AND  TRANET  TRACKING  BIASES  RELATIVE  TO 
A COLLOCATED  LASER 


D.  V.  Carney 
H.  C.  Parker 
J,  H>  Berbert 

As  part  of  the  GEOS-Il  Observation  Systems  Inter  comparison  Investigation, 
the  Wallops  Island  Collocation  Experiment  (WlCE)  was  conducted  during  April  through 
June  1968.  Among  the  collocated  tracking  systems  were  a NASA  laser  and  an 
AN/FPQ-6  C-band  radar  equipped  with  a coherent  signal  processor  (CSP),  The 
FPQ-6  tracked  GEOS  both  in  a beacon  mode  and  in  a skin  tracking  mode.  During 
four  of  the  skin  track  segments,  C-band  doppler  data  was  successfully  taken  and  re- 
corded for  analysis. 

The  laser  data  were  used  to  form  reference  orbits  against  which  the  doppler 
data  were  compared  with  the  NONAME  orbit  dete'*minatlon  program.  Simple  error 
models,  consisting  of  measurement  and  timing  bias  terms,  were  then  fit  to  the  data. 
The  RMS  noise  of  the  residuals  varied  from  2,7  cm/sec  to  13,3  cm/sec,  and  the  de- 
rived range  rate  measurement  bias  varied  from  -1.3  ±0.6  cm/sec  to  6.7  ±0,7  cm/sec 
with  an  average  value  of  2. 5 cm/sec.  The  timing  bias  ranged  from  0,01  ± 0. 20  milli- 
seconds to  0.45  ± 0.20  milliseconds. 
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INTRODUCTION 


The  GEOS-II  Observation  Systems  Intercomparison  Investigation  was  de- 
signed to  evaluate  the  relative  accuracies  of  various  geodetic  observation  systems.  As 
a part  of  this  effort,  the  Wallops  Island  Collocation  Experiment  (WICE)  (Refereoce  1) 
was  conducted  during  the  period  of  April  through  June  1968.  A NASA  laser,  two 
C-band  radars  (an  AN/FPQ-6  and  an  AN/FPS-16),  several  cameras,  a Navy  TRANET 
doppler,  and  an  Army  SECOR  comprised  the  tracking  systems  employed.  These 
trackers  were  collocated  with  one  another  in  order  to  minimize  relative  station  survey 
errors  (determined  to  better  than  lO  cm. ) and  earth  gravity  field  errors  that  might 
otherwise  be  aliased  into  tracking  biases.  The  NASA  laser  range,  azimuth,  and  eleva- 
tion (RAE)  data  were  weighted  at  2 meters,  200  arc  seconds,  and  200  arc  seconds, 
respectively , and  were  used  to  form  the  reference  orbits  against  which  the  other  sys- 
tems were  compared. 

The  AN/FPQ-6  C-band  radar  used  a coherent  signal  processor  (CSP)  to 
provide  doppler  range  rate  measurements  on  GEOS-II  during  the  WICE  (References  2 
and  3).  However,  since  GEOS  did  not  have  a coherent  beacon,  the  CSP  could  only  be 
used  during  a skin  track.  7.'he  skin  tracks  were  augmented  by  a Van  Atta  retro- 
reflector  array  on  the  GEOS  spacecraft.  Of  the  3-1  FPQ-6  tracks  simultaneous  with 
the  laser,  10  were  taken  loth  in  the  beacon  and  skin  track  modes.  For  those  10 
passes,  the  FPQ-6  was  calibrated  for  the  skin  track  portion,  tracked  the  first  third  of 
the  pass  with  the  beacon,  was  switched  to  skin  track  for  the  second  third  of  the  pass, 
and  back  to  the  beacon  for  the  final  third.  Doppler  tracking  was  attempted  during  nine 
passes,  successfully  on  «ight  of  the  nine.  Of  the  eight,  the  four  for  which  we  received 
data  (on  May  30,  June  5,  June  11  and  June  12)  are  covered  in  this  report. 

1.2  SOFTWARE 

The  FPQ-6  doppler  data  forwarded  to  us  fl:om  Wnllops  Station  were  known 
to  contain  a hardware  truncation  error  with  a maximum  value  of  2. 8575  cm/sec.  Sim- 
ilarly, the  zero-set  bias  correction  normally  measured  at  Wallops  had  not  been  applied 
or  recorded  (Reference  4).  Neither  of  these  errors  were  recoverable.  The  data,  as 
received , were  corrected  for  a known  error  in  the  speed  of  light  constant  representa- 
tion which  existed  in  the  radar  hardware  system  at  that  time  (Reference  5)  and  were 
preprocessed  with  the  RCA  C-band  Pre-processing  Program  (Reference  6)  for  time 
tag  and  tropospheric  refraction  corrections. 
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As  in  the  WICE  report  (Reference  1) , the  NONAME  Orbit  Determination 
Program  (ODP)  (Reference  7)  was  used  to  reduce  the  laser  RAE  data  through  a least 
squares  fit  procedure  to  form  the  reference  orbits  against  which  the  FPQ-S  and 
TRANET  range  rate  data  were  compared.  Simple  error  models  for  measurement 
bias  ai^  timing  bias  were  then  fit  to  the  range  rate  residuals  with  the  ODP. 

1,3  A PRIORI  DATA  AND  ERROR  MODEL  WEIGHTINGS 

For  the  purposes  of  this  report,  a specific  jrror  model  is  defined  as  any 
distinct  combination  of  bias  uncertainty  terms  and  the  values  assigned  to  those  terms , 
i.  e. , different  value  combinations  for  the  bias  uncertainty  terms  create  different 
error  models. 

Of  the  four  available  passes , two  were  simultaneous  with  TRANET  doppler 
data,  on  June  11  aol  June  12.  Previous  analyses  of  the  TRANET  data  have  established 
appropriate  weightings  for  that  data  and  confidence  in  the  rebuUs.  It  was  determined 
that,  due  to  the  fairly  high  correlation  of  the  range  rate  bias  to  the  timing  bias , the 
TRANET  time  bias  had  to  be  constrained  to  0.2  milliseconds  (ms)  in  order  to  get 
reasonable  estimates  of  the  measurement  bias.  However,  the  first  FPQ«6  error 
model  allowed  an  increase  to  1 millisecond  in  the  a priori  timing  uncertainty  in  an 
attempt  to  estimate  the  relatively  unknown  timing  bias  for  the  new  C«*band  doppler 
data.  Since  reasonable  estimates  of  the  timing  and  measurement  bias  were  not  re- 
covered, it  was  again  assumed  that  due  to  the  accuracies  of  the  collocation  techniques, 
the  relative  timing  bias  was  as  accurate  pa  that  of  the  TRANET.  Thus , In  the  second 
error  model , the  a priori  FPQ-6  timing  uncertainty  was  reduced  to  0. 2 milliseconds , 
matching  the  TRANET  error  model  terms. 

The  first  two  rows  of  Table  I give  the  a priori  weights  (SIGMA)  assigned  to 
the  laser  measurements  used  to  form  the  reference  orbits.  The  third  row  indicates 
the  SIGMA  applied  to  the  FFQ-G  range  rate  data  in  the  unmodeled  runs  for  each  of  the 
four  passes.  For  both  remaining  stations,  the  columns  show  the  error  model  number, 
the  SIGMA,  the  Initial  estimate  of  the  measurement  bias  (BIAS),  the  uncertainty  about 
that  measurement  bias  (BIAS  SIGMA) , the  initial  estimate  of  the  measurement  timing 
bias  (TIME  BIAS) , and  the  uncertainty  about  that  time  bias  (TIME  SIGMA) . Error 
model  1 was  used  in  the  original  analyses  of  the  FPQ-6  range  ra;e  data;  models  2 and 
3 came  from  the  WICE  TRANET  error  modeling  study. 
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Id  all  of  the  error  model  recovery  runs,  the  laser  reference  orbit  v'as  held 
fixed  by  weighting  the  position  and  velocity  state  vector  components  at  10“®  meters  and 
10“^^  meters  per  second,  respectively,  thus  forcing  the  error  model  coefficients  to 
fit  the  C'band  residuals. 

2.1  RESULTS 

Figure  I is  a display  of  the  plots  of  the  unmodeled  range  rate  residuals 
against  the  laser  RAE  reference  orbits  for  each  of  the  four  passes.  Each  plot  shows 
the  FPQ-6  residuals  sampled  every  five  seconds  and  the  time -location  of  the  doppler 
skin  track  portion  within  the  laser  pass.  None  of  the  FPQ-6  range  rate  initial  points 
are  coincident  with  the  first  range  point  of  the  skin  track;  rather,  the  range  rate  data 
starts  later  by  58,  65,  79  and  109  seconds,  respectively,  for  the  four  passes  taken  in 
chronological  order.  Further,  the  point  of  the  closest  approach  (PCA)  of  the  satellite 
consistently  occurs  at  the  second  time  nolnt  of  each  range  rate  track  and  lies  between 
71®  and  88®  elevaUon. 

The  lower  two  plots,  i.e. , the  June  11  and  12  passes,  also  show  the  un- 
modeled TRANET  residuals  for  comparison.  These  data  were  smoothed  over  a 32- 
second  interval  during  d;e  Naval  Weapons  Laboratory  (NWL)  preprocessing  and  are 
shown  without  further  sampling. 

Figure  II  illustrates  the  residuals  as  they  appear  after  the  second  error 
model  was  applied  during  t h the  FPQ-6  and  TRANET  data  reductions. 

Table  n summarizes  the  pass  statistics  for  the  residuals  and  the  derived 
measurement  and  timing  bias  error  model  terms  for  the  four  FPQ-6  tracks  and  the 
two  simultaneous  TRANET  tracks.  The  RMS  noise  values  are  given  for  both  the  first 
ana  final  iterations  to  show  the  effective  noise  reduction  resulting  from  modeling  the 
systematic  trends  within  each  pass.  The  mean  value  of  the  residuals  for  iteration  one 
is  given  for  comparison  to  the  measurement  bias  term  derived  in  subsequent  iterations, 
when  the  residual  mean  has  been  reduced  to  zero.  When  both  the  biases  and  their 
associated  uncertainties  are  considered,  the  best  results  appear  to  be  obtained  with 
the  second  error  model. 

The  third  error  model  was  used  in  an  attempt  to  display  the  effects  caused 
by  assuming  the  measurement  bias  to  be  fairly  well  known  and  the  timing  bias  relatively 
unknown. 
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Table  HI  is  a summary  of  the  statistics  for  each  item  of  interest  averaged 
over  the  four  FPQ-6  passes,  the  two  simultaneous  passes  of  the  low  frequency  TRANET 
pair  <Lo) , and  finally,  for  reference,  over  the  26  passes  of  the  low  and  the  16  passes 
of  the  high-frequency  TRANET  pairs  taken  during  the  WICE, 

It  has  been  shown  in  a separate  set  of  runs  not  given  here , that  when  no 
time  bias  error  term  is  solved  for,  the  derived  measurement  bias  for  the  final  iteration 
is  equal  to  the  residual  mean  value  from  the  first  iteration.  The  measurement  errors 
are  forced  into  the  range  rate  bias  with  a slight  degradation  (0.0  to  0. 1 cm/sec.)  in 
the  RMS  noise;  i.e. , if  a timing  bias  does  indeed  exist,  it  can  be  aliased  into  a meas- 
urement bias  with  an  appropriate  error  model.  Conversely , if  no  measurement  bias 
term  is  allowed  and  the  time  bias  uncertainty  is  set  large  at  ± 1 second,  the  residual 
mean  values  are  reduced  to  a near-zero  level.  The  time  bias  term  will  now  absorb 
most  of  the  measurement  error,  varying  from  -0.30  milliseconds  to  2.51  milliseconds. 
This  indicates  that  the  range  rate  and  timing  biases  are  highly  correlated  and  points 
out  the  need  to  assign  relative  uncertainties  w'hich  approximate  reality.  Due  to  the 
careful  relative  time  synchronization  techniques,  we  believe  that  the  clocks  were 
indeed  synchronized  to  at  least  200/xsecs.  (0,2  millisec.),  as  used  in  error  model  2. 

The  C-band  doppler  tracks  are  all  one-sided,  starting  a few  seconds  before 
PCA  and  then  continuing  down  to  an  elevation  of  between  45°  and  60°.  The  doppler 
data  were  offset  from  the  start  of  the  skin  track  due  to  difficulties  in  "locking  on"  for 
range  rate.  Had  the  passes  been  more  nearly  symmetrical,  the  bias  terms  should 
have  been  more  separable. 

Figure  III  consists  of  two  overlays  of  the  FPQ-6  range  rate  biases  on  plots 
of  the  WICE  TRANET  measurement  biases.  In  both  cases,  the  TRANET  range  rate 
biases  were  derived  with  the  second  error  model.  The  first  (top)  plot  shows  the 
FPQ-6  biases  as  they  were  originally  derived  with  the  first  error  model  (Reference  8), 
The  second  (bottom)  plot  compares  the  FPQ-6  biases  as  derived  with  the  second 
model  to  the  TRANET  biases.  It  is  apparent  that  the  more  realistic  second  model 
yields  FPQ-6  measurement  biases  that  generally  agree  more  closely  with  the  laser 
reference  orbit  and  the  TRANET  than  does  the  first  model. 

Figure  IV  is  an  overlay  of  the  FPQ-6  skin  track  range  biases  on  plots  of 
the  FPQ-6  beacon  mode,  the  FPS-16  C-band,  and  the  SECOR  range  biases  from  the 
WICE,  It  is  Included  as  supplementary  information,  showing  the  degree  of  aggreement 
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H between  the  ranging  biases  for  the  skin  tr'ick  portion  (X)  and  the  beacon  track  portion 

(V ) with  the  laser  reference  orbit  and  the  other  systems  during  that  period. 

fl  U has  been  reprrted  (Reference  9)  that  the  phase  center  of  the  TRANET 

* antenna  was  on  the  order  of  1 meter  different  from  the  geometric  antenna  center  that 

^ we  have  been  using.  Runs  were  made  for  the  June  11  and  June  12  TRANET  passes 

I using  the  modified  phase  center  coordinates.  The  measurement  biases  recovered 

using  error  model  2 deviated  from  the  previous  results  by  less  than  0.4  cm/sec. 

I The  timing  biases  did  not  change. 
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CONCL  JSIQNS 


While  it  is  understood  that.,  due  to  the  sparse  number  of  passes,  these  re- 
sults are  not  very  conclusive,  they  are  an  Indication  of  what  was  achieved  in  early 
GEOS -II  attempts  to  utilize  CSP  doppler  data  availab.le  from  the  Wallops  AN/FPQ-6 
radar  operating  in  a skin  track  mode.  The  data  in  this  study  is  somewhat  corrupted 
with  non-recorded  zero-set  and  truncation  errors , yet  the  results  appear  reasonable 
for  a C-band  range  rate  system  at  that  time,  particularly  for  the  last  (6/12)  pass. 

Considering  error  model  2 , the  FPQ-6  range  rate  biases  range  from 
-1.3  cm/sec  to  6. 7 cm/sec  with  timing  biases  of  0,01  milliseconds  to  0,45  milli- 
seconds . For  the  two  simultaneous  TRANET  passes , the  range  rate  biases  are 
1.9  cm/sec  and  1,2  cm/sec  with  timing  biases  of  0,00  ms.  and  0,05  ms. 

The  measure  me  rt  range  rate  bias  results  are  significantly  better  with 
less  uncertainty  when  the  time  bias  term  is  more  tigntly  constrained  (from  1 second, 
or  1 millisecond,  down  to  0,2  millisecond  uncertainty) , indicating  that,  due  to  the 
correlation  betw'een  the  range  rv\te  and  timing  bias  coefficients,  it  Is  important  to  use 
realistic  time  bias  constraints.  By  using  the  more  realistic  estimate  of  the  uncertainty 
in  the  FPQ-6  clock,  the  results  are  substantially  improved  over  previous  analyses. 
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AVERAGED  STATISTICS  FOR  THE  RANGE  RATE  PASSES 


Sample  RMS  ^ 

Si.e  (cm/sec) 


10. 5 ± 4. 4 
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FIGURE  32.  WICE  RANGE  BIASES  VS.  DATE 
(LASER  REFERENCE  ORBITS) 
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COMPANION  DATA  SETS  WITHOUT  P WERE  PROCESSED 


FLIGHT  TIME  (Seconds) 


FLIGHT  TIME  (Seconds) 


+ 

+ 

* 

+ 

LU 

+ 

cr 

+ 

a 

3 

+ 

o 

LU 

to 

-J 

+ 

+ 

UJ 

=5 

< 

+ 

to 

Ll. 

H 

< 

+ 

H- 

<c 

<C 

o 

K- 

Ct 

< 

a. 

O 

*0;: 

LU 

UJ 

< 

< 

>• 

K 

az 

tx 

o 

to 

CO 

t— 

LU 

LU 

o 

■f 

h- 

h- 

LU 

+ 

O 

O 

+ 

z 

z 

<c 

+ 

LU 

LU 

QC 

+ 

O 

O 

h- 

+ 

+ 

+ 

• 

(saiirt  ivoiinvN)  aoaaa  lovdwi 


FLIGHT  TIME  (Seconds) 


oc 

C3  =, 

o c/»  — 

uj  r>  < 

- li. 

< < o 

t—  o ;r 

Q .(C  ^ 

LU  UJ  ' 


tC  Ofi  “■ 

C/>  C/)  H- 

LU  LlJ  O 

H-  I — 1 1 1 

o o -> 

z z <c 

uj  lu  oe 

Q Q »— 


(S31IW  ivoiinvN)  yoaa3  iovdwi 


■V*-»  » ;* 


FLIGHT  TIME  (Seconds) 
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RELATIVE  TO  A GIVEN  FLIGHT  TIME,  R,A,E,R  DATA  FROM  3 SITES 
WERE  PROCESSED  IN  THE  ORDER  SHOWN  BELOW. 
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SVAFB  TPQ-18  R Versus  INERTIAL  GUIDANCE.  OPERATION  3782,  11  June  1971 
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DOPPLER  COVERAGE,  MINUTEMAN  111  OPERATIONS 
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DOPPLER  COVERAGE,  MINUTEMAN  111  OPERATIONS 


INERTIAL  GUIDANCE  AND  RADAR  TIMING  ERROR 


AT  SVAFB  TPQ-18  and  TPQ-18  R Versus  INERTIAL  GUIDANCE. 
OPERATIONS  3782  and  6448 


POINT  MUGU  FPS-16^4  R Venus  INERTIAL  GUIDANCE.  OPERATION  5411 


SAN  NICOLAS  ISLAND  FPS-!6*3  R Versus  INERTIAL  GUIDANCE,  OPERATION  5179 
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COHERENT  TRACKING  OF  MINUTEMAN  III 
MEASUREMENT  SYSTEM  PROBLEMS  AND  EVALUATION 


WILLIAM  COLLINS.  FEC 

Federal  Electric  Corporation  ITT ' 
Vandegberg  AFB,  California  93437 


SAN  NICOLAS  013003  FPS-16  VESS  INSTALLED  OCTOBER  1968 
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IMPROVED  DATA  TO  RANGE  USERS 


LOW  BEACON  TRANSMITTER  POWER 
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G40  PRF  REDUCES  EFFECT 


65  Seconds  SIGNAL  SPECTRUM 
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*3.0 


ROOT  Vs  IG  (BET  AUTO) 


FLIGHT  TIME  (Seconds) 
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MINUTEMAN  III  EVALUATION 
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PROVIDES  INCREASED  CONFIDENCE  IN  EVALUATION  RESULTS 
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SUMMARY  OF  TEST  RESULTS 
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PVM-3  8/24/73  N/A  ~1  fps  POWERED 
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PROBLEMS  ASSOCIATED  WITH  DEVELOPMENT  OF  A LAUNCH  HEAD  RANGE  SAFETY 
SYSTEM  FOR  CONTAINING  ICBM  LAUNCHES  IN  THE  KWAJALEIN  LAGOON  CORRIDOR 

Stan  Radom,  SAMTEC 

{Verbatim  Transcript  from  Tape) 


I am  using  a viewgraph  made  in  1966  in  order  to  give  you  an  idea  of  vhy  we  got  into 
coherent  tracking  in  the  first  place.  There's  some  interesting  history  attached  to  it.  I'd  like 
to  go  over  briefly  why  the  very  large  expenditures  were  made  here,  San  Nicolas  Island, 
Pillar  Point,  and  Point  Mugu  to  get  into  the  range  and  range  rate  business. 

When  the  Western  Test  Range  was  formed,  the  Air  Force  took  over  PMR  responsibilities 
for  Instrumentation  of  ballistic  launches  and  space  launches  from  the  west  coast.  The 
Electronics  System  Division  of  the  Air  Force  was  assigned  the  task  of  attempting  to 
outguess  our  needs.  They  were  going  to  tell  us  what  we  were  going  to  need  to  accomplish 
our  objectives. 

Back  in  about  1966  there  was  a detachmen.  of  ESD  people  here.  When  they  took  a 
look  at  this  problem  they  said,  "Yes,  you  are  suffering  from  a lack  of  downrange  islands  so 
something  drastic  has  to  be  done."  The  requirement  was  thrust  upon  us  to  start  impacting 
ICBMs  in  the  Kwajalein  area.  That  was  because  a great  deal  of  money  was  being  spent  in 
the  Roi'Narnur  complex  for  TRADEX  and  other  radars  associated  with  evaluation  of 
reentry  systems.  Then  they  added  the  requirement  to  score  impacts,  and  recover  warheads. 
The  recovery  dictated  the  need  for  the  lagoon  impacts.  They  not  only  wanted  an  accurate 
score,  they  also  wanted  to  get  the  warhead  (or  pieces)  back.  That  dictated  negotiation  with 
the  Marshallese  to  evacuate  a corridor  32  miles  wide  through  the  middle  of  the  Kwajalein 
lagoon. 


(See  Viewgraph) 

This  basic  problem  and  why  I brought  it  up,  points  out  the  need  and  justifies  the 
expenditure  for  the  coherent  modifications  to  the  FPQ-6  at  Pillar  Point  and  the  FPS-16's  at 
PMR.  The  TPQ‘18  was  already  programmed  for  coherent  mod  by  the  ETR  before  it  was 
delivered  to  Vandenberg. 

So  this  is  generally  the  problem.  A much  different  situation  because  you  have  natives, 
plus  military  and  civilian  dependents  on  Kwajalein  and  you  have  a population  up  at 
Roi'Narnur.  This  was  drastically  different  than  the  situat.:in  at  Eniwetok  where  there  were 
no  natives,  no  dependents,  and  where  shelter  was  provided  for  all  essential  personnel.  It  was 
a good  place  to  do  extra-hazardous  testing.  At  that  time  we  had  developed  an  underwater 
scoring  system  for  the  Eniwetok  lagoon  that  exceeded  our  expectations,  and  we  had  a 1 
Sigma  circular  accuracy  of  34  feet  on  impacts.  Of  course,  you  can't  argue  when  the 
recovery  barge  is  vectored  into  position,  and  the  diver  dives  down  about  200  feet  and  steps 
on  the  warhead  I But  with  this  requirement  at  Kwajalein,  we  had  to  do  something  drastic, 
so,  ESD  initiated  a study  with  the  help  of  MITRE,  to  see  if  multilateration  radars  bearing 
down  on  the  trajectory  might  provide  a range  safety  system  that  could  contain  the  impacts 
in  the  corridor  without  inadvertent  destruction  of  good  missiles. 
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So  we  took  a look  at  what  the  geometrical  configuration  could  do.  We  had  some 
simulation  software;  we  had  a contract  with  RCA  to  take  a look  at  this  situation,  along 
with  the  work  being  done  by  ESD,  and  what  little  work  we  did  in-house  at  the  time.  This 
generally  is  wnat  the  situation  looks  like. 


{First  Overlay) 

With  a single  FPS-16  at  the  launch  head,  and  a computer,  you  have  an  uncertainty  of 
something  like  42  miles,  bigger  than  the  corridor  itself  — its  just  out  of  the  question!  The 
future  impact  prediction  was  just  not  useful.  There  was  no  question  but  that  it  was 
inadequate.  The  multilateration  mode,  that  I'm  going  to  show  you  here,  closely  correlates 
both  the  ESD  study  and  RCA  analysis,  and  what  everyone  took  a look  at. 

i.ets  take  a look  at  adding  radars,  and,  although  we  don't  have  downrange  islands,  lets 
: /lat  wc  can  do  with  bilateration. 


jr.:oxt  Overlay) 

This  reduced  the  IIP  uncertainty  in  ordar  of  something  like  10  miles.  Still 
unsatisfactory  • it  didn't  allow  enough  latitude  for  the  excursion  that  a MINUTEMAN,  for 
example,  might  make  in  the  last  few  seconds  of  powered  flight.  Adding  another  sensor  — If 
you  had  conventional  pulse  radars.  Point  Pillar  and  San  Nicolas  Island  giving  you  the  longest 
baseline  for  bilateration  plus  trilateration,  using  a ship  providing  a data  link,  as  you  can  see 
you  then  have  an  uncertainty  el  ipse  that  looks  tike  your  improving  things.  Of  course:,  the 
ship  can  give  you  the  Z-component,  measuring  altitude  in  an  optimum  position  near 
burnout.  Of  course,  its  the  Z-component  when  you  go  a quarter  of  the  way  around  the 
earth,  4300  mites  to  the  target  area  at  Kwajatein,  thats  your  major  source  of  error.  So  you 
can  see  why  there  was  a dramatic  improvement  here,  but  it  was  still  not  satisfactory.  The 
risk  would  be  too  great  to  say  that,  if  we  had  conventional  pulse  radars,  bilateration,  plus  a 
ship  and  real  time.  So  then  we  took  a look  at  what  range  and  range  rate  would  do  for  you. 


(Next  Overlay) 

As  soon  as  you  went  coherent  in  a bilateration  mode,  when  reading  range  and  range 
rate  directly,  you  can  see  that  the  size  of  the  uncertainty  ellipse  narrowed  down,  moved 
about  90  degrees.  But  when  /ou're  bouncing  back  and  forth  against  the  corridor  lines,  it 
still  hasn't  done  enough  good. 

Lets  take  a look  at  what  happens  when  you  put  a ship  in  the  solution.  This  is  a ship 
near  burnout  with  a conventional  pulse  radar,  but  providing  real-time  data  reliably  into  th(3 
computation  center.  It  reduces  your  uncertainty  ellipse  to  something  like  2 by  3 miles.  And 
that  boks  goodi  Theory  showed  that  this  is  the  way  to  gol 

Lets  go  coherent  at  Point  Pillar  and  give  the  Navy  i.he  money  to  modify  radars  at  San 
Nick,  Point  Mugu,  and  awuy  we  gol 

In  those  days  we  had  ships  also.  This  convinced  the  people  reviewing  the  budget,  and 
everyone  else,  that  this  was  a sound  method  of  approach.  All  the  problems  you've  been 
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hearing  about,  some  of  the  presentations  that  our  folks  have  shown  you,  kept  talking  about 
the  I IPs  and  Kwajalein,  this  was  the  compelling  requirement.  The  problem  has  been  taking 
the  theory  and  reducing  it  to  practice. 

As  in  MINUTEMAN,  the  problem  of  having  a range  safety  system  for  a multi-stage  solid 
rocket,  with  high  accelerations  is  that  it  requires  that  action  be  taken  in  the  last  half-second 
of  powered  flight.  You  have  to  measure  lateral  acceleration  near  burnout  to  reliably  contain 
the  impact  within  the  boundaries  of  Kwajalein  lagoon  corridor.  That  means  you  are 
going  to  have  to  have  computer  abort.  You  can't  have  a guy  at  a console  with  his  finger  on 
a button.  The  brain-to-finger  response  time  alone  is  out  of  the  question. 

As  soon  as  we  started  talking  about  computer  abort,  the  Range  Commander  who  has 
the  safety  responsibility  and  can't  delegate  it,  asked  the  great  question:  "How  reliable  is  my 
system,  what  is  the  possibility  of  that  computer  blowing  a good  'bird'?" 

So,  my  small  analysis  group  started  looking  at  reliability  figures  for  radars  in  general, 
data  links,  everything  associated  with  getting  the  data  into  the  computer.  Assuming  that  we 
had  the  software  that  worked  at  the  time,  what  were  the  chances  of  a dropout  that  would 
cause  a computer  abort?  It  turns  out  that  the  Range  Safety  Officer  would  buy  about  5 
percent  of  the  birds.  So,  each  Commander,  when  given  that  information,  said  "I'm  going  to 
take  my  chances,  I don't  want  computer  abort.  I'd  rather  go  with  the  waiver  before  I design 
into  the  system  an  assured  method  of  blowing  1 out  of  20  'birds'."  Generals  Bleymaier, 
Kronauer,  Wilson  and  Lowe  felt  the  same  way.  We  haven't  had  a Commander  who  has  felt 
differently.  As  a result,  this  kind  of  put  the  haitus  on  where  we  were  going  with  this 
system.  Furthermore,  we  didn't  have  the  computational  power,  as  an  example,  even  If  we 
had  a system  that  could  work  reliably.  Because  we  did  not  demonstrate  a positive  system  of 
range  safety,  we  couldn't  tell  our  users,  "you  must  carry  a coherent  beacon." 

That  is  essentially  where  we  are  today,  and  one  of  the  reasons  for  this  conference.  We 
have  coherent  mods.  We  know  a lot  more  about  beaconry  and  we  know  a lot  more  about 
the  reliability  of  data  links.  So  we  need  another  appraisal  - .maybe  an  agonizing  reappraisal 
of  whether  we're  going  to  take  another  stab  at  this  system,  take  a look  at  GERTS  for 
MINUTEMAN,  or  whatever  means  we  can.  Maybe  arj  interferometric  system?  I'm  not 
talking  about  a floating  MISTRAM  or  something  like  that,  but  by  taking  another  look  at 
whatever  technology  is  available  to  us,  so  we  can  go  from  the  waiver  mode  to  a positive 
system  of  range  safety  for  Kwajalein. 

I hope  that  this  briefing  will  give  you  some  background  on  why  the  initial  investment 
was  made  in  coherent  radars  for  the  West  Coast  launch  head. 
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PULSE  DOPPLER  INSTRUMENTATION  RADARS 
Renzo  Mitchell,  RCA 

FPQ-6/TPQ-18 

VAFB 

PILLAR  POINT 
WALLOPS  ISLAND 

FPS-16 

PT.  MUGU 
SAN  NICOLAS 

MPS-36 

2- SAMTEC  (MRSS) 

1 - TONOPAH 

2 - KWAJALEIN 

3 - GERMANY 


6 - WSMR 


DESIGN  PERFORMANCE  CAPABILITY 


FPQ-6/TPQ-18 

RANGE  RATE  0.1  FT./SEC. 

COHERENT  INTEGRATION  IN  RANGE  AND  ANGLES 

FPS-16 

RANGE  RATE  0.1  FT./SEC. 

MPS-36  1.0  FT./SEC.  (SPEC.) 

0.1  FT./SEC.  DESIGN  CAPABILITY 

MPS-36  FINE  LINE  TRACK 

FINE  LINE  TRACK  IN  RANGE  AND  ANGLES. 

DECREASED  DOPPLER  AMBIGUITY  RESOLUTION  TIME. 


CAPABILITY  OF  TRACKING  THROUGH  CLUTTER. 


IMPROVEMENTS  IN  CAPABILITY 


SOFTWARE  ENCHANCEMENT  TO  MAINTAIN  DOPPLER  TRACK  DURING  STAGING. 

SOFTWARE  ENHANCEMENT  TO  MAINTAIN  DOPPLER  TRACK  DURING  EVENTS  THAT 
CAUSE  SPECTRAL  SPREADING  - FLAME  EFFECTS. 

SATELLITE  COMPUTER  TO  RELIEVE  4101  OR  DDP-124. 

AUTO  CHECKOUT  TO  ASSURE  OPERATIONAL  READINESS. 

VELOCITY/ACCELERATION  PROCESSING  IMPROVEMENTS  TO  ENHANCE  INITIAL 
LOCK-ON  CAPABILITY  AT  LOWPRF  AND  HIGH  ACCELERATION. 

DRIVER  STROBE  RISE-TIME  REDESIGN. 

SYSTEM  NOISE  FLOOR  IMPROVEMENTS. 
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GENERAL  PROJECT  DESCRIPTION 


The  purpose  of  the  GEOS-C  Project  is  to  design,  develop,  and  launch 
a geodetic  and  oceanographic  satellite  and  to  perform  experiments  in  support 
of  the  application  o::  geodetic  satellite  techniques  to  the  Earth  and  Ocean 
Physics  Applications  Program, 

The  GEOS-C  Mission  will  provide  data  with  which  to  refine  the  geodetic 
and  geophysical  results  of  the  National  Geodetic  Satellite  Program  (NGSP) 
and  will  furnish  a test  bed  for  new  systems.  This  Mission  will  also 
contribute  to  fulfilling  the  C-band  radar  calibration  and  altimeter  re- 
quirements of  the  Departments  of  Defense  and  Commerce. 

The  spacecraft  for  this  mission  was  designed  and  fabricated  by 
the  Applied  Physics  Laboratory  of  the  Johns  Hopkins  University.  The 
structural  configuration  is  based  on  the  GEOS-2  mechanical  design. 

Basically,  the  structure  will  be  the  same  as  GKOS-2  with  the  sub- 
stitution cf  heavier  trusses  to  accommodate  the  additional  weight. 

The  GEOS-C  (Geodynamics  Experimental  Ocean  Satellite)  will  be 
launched  during  FY-75  from  the  Air  Force  Western  Teat  Range  (AFWTR) 
located  at  Vandenberg  AFB,  California.  TChe  nominal  orbit  parameters  are; 
Mean  Altitude  - 843  km 
Inclination  - 115  degrees 
Eccentricity  - 0.006 

- 101.8  minutes 
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Period 


The  orbital  parameters  are  chosen  to  maximize  the  possibility  of  pro- 
viding orbit  traces  which  cover  the  earth’s  surface  in  a prescribed 
grid  work  pattern, 1. e. , 1®  x 1®. 

Following  orbital  injection,  the  spacecraft  will  employ  an  attitude 
stabilization  system  consisting  of  a gravity  gradient  boom  with  end  mass 
and  damper,  momentum  wheel  and  an  electromagnet.  The  energized  electro- 
magnet creates  a magnetic  dipole  moment  along  the  boom  axis  of  the 
satellite  to  aid  in  gravity  gradient  capture.  The  axis  will  align  with 
the  earth's  magnetic  field  so  that  in  the  northern  latitudes  the  proper 
side  of  the  spacecraft  will  point  earthward.  At  this  time,  the  boom  will 
be  extended  and  the  electromagnet  turned  off. 

Power  for  operation  of  the  spacecraft  and  experiment  sub-systems 
will  be  provided  by  a single  11-cell  battery.  The  battery  will  provide 
a 14.7  volt  nominal  power  supply  for  the  duration  of  the  mission  which 
is  scheduled  for  a period  of  one.  year. 

A PCM/FM  telemetry  system  will  allow  two  basic  modes:  low  bit 

rate  data  transmission  (1.562K  bits/sec)  or  high  bit  rate  data  transmission 
(15.62K  bits/sec.).  Transmission  of  data  will  be  either  via  S-band  direct 
to  ground,  S-band  to  ground  via  the  ATS-F  satellite,  and/or  VHF  direct  to 
ground.  Experiment  sub -systems  will  be  controlled  by  a series  of  commands 
either  stored  (delayed  command)  or  sent  directly  from  ground  stations. 


GEOS-C  MISSION  OBJECTIVES 


The  GEOS-C  Mission  Obj actives  In  order  of  priority  at  launch 

are: 


To  perform  an  in-orbit  satellite  altimeter  experiment  to;  (1) 
determine  the  feasibility  and  utility  of  a space-borne  radar  altimeter 
to  map  the  topography  of  the  ocean  surface  with  an  absolute  accuracy 
of  meters,  and  with  a relative  accuracy  of  1-2  meters,  (2)  determine 
the  feasibility  of  measuring  the  deflection  of  the  vertical  at  sea, 

(3)  determine  the  feasibility  of  measuring  wave  height,  and  (4)  contribute 
to  the  technology  leading  to  a future  oner  ^tional  altimeter-satellite 
system  with  a 10-cm  measurement  capabil; 

* To  further  support  the  calibration  of  NASA  and  other  agencies* 
ground  C-band  radar  systems  by  providing  a space-bome  'oherent  C-Band 
transponder  system,  to  assist  in  locating  these  stations  in  the  unified 
earth- centered  reference  system,  and  to  provide  tracking  coverage  in 
support  of  the  radar-altimeter  experiment. 

• To  perform  a satellite-to-satelllte  experiment  Ith  the  ATS-F 
satellite  using  an  S-Band  transponder  subsystem  to  dir^'ctly  measure  the 
short  period  accelerat:'  imparted  to  the  spacecraft  by  the  gravity 
field  and  to  determine  the  position  of  the  spacecraft.  The  anticipated 
measurement  data  quality  of  about  .04  cm/sec  over  a ten-second  integration 
interval  will  aid  in  Improving  the  earth  gravity  model  up  to  spherical 
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harmonic  terms  of  degree  and  order  to  approximately  25  and  in  providing 
tracking  coverage  over  mid-ocean  areas  to  support  the  radar-altimeter 
experiment*  The  satellite-to-satellite  system  will  also  be  used  for 
altimeter-data  relay  through  AFS-F. 

* To  further  support  the  inter comp aria on  of  new  and  established 
geodetic  and  geophysical  measuring  systems  including:  the  radar  altimeter, 

satellite-to-satellite j c-Band,  S-Band,  laser,  and  doppler 
tracking. 

To  investigate  solid-earth  dynamic  phenomena  such  as  polar 
motion,  fault  motion,  earth  rotation,  earth  tides,  and  continental  drift 
theory  with  precision  satellite  tracking  systems  such  as  laser  and 
doppler  ground  stations* 

To  further  refine  orbit-determination  techniques,  the  determination 
of  interdatum  ties,  and  gravity  models  with  a spacecraft  equipped  with 
laser  retroref lectors,  C-Band  and  S-Band  transponders  and  doppler  beacons. 

To  support  the  calibration  of  S-Band  sites  in  the  STDN  by 
furnishing  a space-borne  transponder,  to  assist  in  positioning  the 
network  stations  in  the  world  reference  system,  and  to  assist  in  evaluating 
the  system  as  a tool  for  geodesy  and  precision  orbit  determination. 

The  achievement  of  these  objectives  not  only  will  constitute  a 
successful  mission  but  also  will  greatly  enhance  man’s  understanding  of 
the  physical  properties  of  the  plar.  -t  Earth, 
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EXPERIMENT  PACKAGE 


The  GEOS-C  experlrsnc  package  will  consist  of  five  basic  instruments 
each  of  which  will  contribute  to  one  or  more  of  the  mission  objectives. 
Consistent  with  these  objectives,  the  experiment  sub-systems  are  listed 
in  priority  order  as  follows; 

Radar  Altimeter 

Coherent  C-Band  Transponder 

S-Band  Instrumentation  for  Satelllte-to-Satelllte  Experiments 
Laser  Retroreflector 

Doppler  Transmitters 
Non-Coherent  C-Band  Transponder 

S-Band  lnstrumentatlo.i  for  Earth  Tracking  Experiments 


Radar  Altimeter 

The  basic  objectives  of  the  radar  altimeter  experiment  are  to  demonstrate 
the  feasibility  of  utilizing  an  on-board  altimeter  to  measure  the  time- 
varying  behavior  of  the  ocean's  surface  and  the  departure  of  the  sea  surface 
from  the  geold  and  to  Investigate  altimeter  Instrumentation  technology.  To 
meet  these  objectives,  the  altimeter  will  have  two  distinct  data  gathering 
modes:  Long  Pulse  and  Short  Pulse.  The  basic  measurement  goals  established 

for  the  altimeter  are  as  follows: 

Precision:  Short  Pulse  Mode  30  cm 

Long  Pulse  Mode  60  cm 

Geold  Accuracy:  Absolute  + 5 meters 

Relative  + 2 meters 

Sea  State:  25%  of  S.V.H. 

One  of  the  primary  goals  of  the  experiment  will  be  to  obtain  engineer- 
ing data  on  altimeter  performance.  For  this,  It  will  be  necessary  to  measure 
and  evaluate  parameters  such  as  sea-surface  roughness  and  spacecraft 
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libration  as  it  relates  to  hardware  performance.  Another  goal  will  be 
to  calibrate  the  altimeter  over  an  ocean  area.  The  presently  planned 
calibration  area  will  be  the  portion  of  the  Atlantic  Ocean  bounded 
by  Wallops  Island,  Virginia;  Merritt  Island,  Florida;  Grand  Turk;  and 
Bermuda.  Altimeter  accuracy  will  be  determined  by  comparing  the 
altitude  measured  by  the  altimeter  to  the  spacecraft  altitude  determined 
by  independent  tracking  systems.  Precision  and  resolution  will  be 
determined  by  comparing  sea  surface  profiles  resulting  from  altimeter 
measurements  to  profiles  determined  by  independent  methods, 

C~Band  System 

Two  C-Band  radar  transponders  will  be  flown  on  the  GEOS-C  satellite 
to  support  the  altimeter  and  C-Band  system  calibration  as  well  as 
geometric,  gravimetric,  and  other  geodetic  investigations.  The  C-Band 
System  consists  of  the  two  transponders  (one  coherent  and  one  non- 
coherent) and  the  associated  ground  tracking  C-Band  radars. 

The  non-coherent  transponder,  operating  in  conjunction  with  existing 
ground  radar  systems,  will  provide  for  range  and  angle  measurements. 

Tlie  coherent  transponder,  operating  in  conjunction  with  existing  coherent 
ground  radar  systems,  will  provide  for  range,  range  rate,  and  angle 
measurements. 

Laser  System 

The  Laser  System  consists  of  the  GEOS-C  spaceborne  laser  retroref lector 
subsystem  and  the  ground-based  Laser  Ranging  Systems.  The  retroreflector 
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will  be  utilized  in  conjunction  with  the  ground-based  laser  systems 
to  obtain  precision  satellite  ranging  data. 

In  the  GEOS-C  time  frame,  the  NASA  Laser  Ranging  Systems  are 
expected  to  have  a ranging  capability  of  10  cm  or  better.  The  capability 
of  the  angle  data  is  estimated  to  be  about  0.5  milllradian  or  better. 

Range  and  angle  data  are  provided  at  a once-per- second  rate. 

Doppler  System 

The  Doppler  System  consists  of  two  spacebome  transmitters  and 
ground  doppler  receiving  stations.  The  dual  frequencies  (162  and  324  MHz) 
originate  from  an  ultra-stable  5 MHz  oscillator* 

Ground  observation  stations  measure  the  doppler  components  of  the 
signals  received  from  the  GEOS-C  spacecraft  by  counting  cycles  resulting 
from  the  difference  between  the  received  frequency  and  the  station 
oscillator.  The  difference  frequencies  between  the  higher  and  lower 
received  frequencies  and  the  station  oscillator  are  combined  in  tne 
proper  proportions  to  obtain  both  the  first  order  ionospheric  refraction 
correction  and  the  refraction  corrected  doppler  frequency. 

S-Band  System 

The  S-Band  system  consists  of  the  following  elements: 

GEOS-C  Coherent  S-Band  Transponder 
GEOS-C  S-Band  Antenna  System 
ATS-F  Spacecraft 
ATS-F  Ground  Terminals 


STDN  S-Band  Ground  Stations 


MISSION  PROFILE 


The  GEOS-C  Mission  has  been  divided  into  two  distinct  phases. 

Phase  I covers  all  activities  from  launch  through  one  year  of  Experiment 
data  collection  and  Phase  II  covers  those  activities  after  Phase  I 
through  the  remainder  of  the  Mission  lifetime. 

Phase  I can  be  sub-divided  into  several  phases  according  to  the 
extent  of  experiment  data  collection , the  type  of  data  being  collected » 
and  various  other  operational  and  physical  constraints.  These  sub-phases 
along  with  the  dominant  activity  are: 


PHASE  1 PERIODS 

TIME  PERIOD 

DOMINANT  ACTIVITY 

(Days  After  Launch) 

Phase  A 

0-40 

Launch  and  Operational 
Assessment 

Phase  B 

40  - 130 

Experiment  Systems  Cali- 
bration & Evaluation 

Phase  C 

130  - .'•85* 

Unique  Experiments  and 
Localized  Grid  Activity 

Phase  D 

285*-  325* 

Global  Activities 

Phase  E 

325*-  405 

Localized  Grid  Densification 

* The  length  of  Phase  D Is  not  expected  to  change.  However,  the  time 
of  occurrence  is  dependent  upon  the  time  of  ATS-F  drift. 


PHASE  A " Launch  and  Operations  Assessment 


This  phase  begins  with  Launch  and  extends  over  a period  of  approx- 
imately 40  days  in  which  the  following  activities  occur: 

(a)  Launch 

(b)  Orbit  Injection 

(c)  Early  orbit  determination  and  refinement 

(d)  Gravity  gradient  capture  and  damping 

(e)  Momentum  wheel  turn  on 

(f)  Yaw  capture  and  stabilization  damping 

(g)  Spacecraft  functional  and  electrical  checkout 

(h)  Operational  assessment  of  experiment  systems 

It  is  not  expected  that  any  useful  experiment  data  collection  will 
be  accomplished  during  Phase  A. 

PHASE  B - Experiment  Systems  Calibration  and  Evaluation 

It  is  expected  that  this  phase  will  begin  about  40  days  after  launch 
and  continue  for  two  or  three  months.  It  is  expected  that  the  bulk 
of  the  data  collected  during  this  period  will  be  useful  for  investigation 
purposes.  However | data  distribution  during  this  period  will  be  slower 
than  normal  due  to  the  more  detailed  analyses  required  to  calibrate  all 
experiment  systems  and  tc  validate  all  data  processes. 

Major  activities  associated  with  this  phase  include  the  following: 

(a)  Altimeter  experiment  systems  calibration 

(b)  Satellite-'to-satelllte  experiment  systems  calibration 
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(c)  Altimeter  experiment  a/stems  operation  globally, within  the 
real-time  ground  TM  station  covera^'^e  areas,  on  those  days  not 
utilized  for  Altimeter  Calibration  activities, 

(d)  Ground  tracking  system  (Laser,  C-Band  and  Doppler)  data 
collection  activities  on  a global  basis  to  the  maximum  extent  possible 
commensurate  with  power  budget » other  systems  calibration  activities,  and 
Investigator  needs. 

PHASE  C - Unique  Experiments  and  Localized  Grid  Activities 

Phase  C activities  are  expected  to  begin  approximately  130  days  after 
GEOS-C  launch  and  continue  until  AXS-F  is  maneuvered  from  the  Western 
Hemisphere  location  (94  degrees  west  longitude)  to  the  Eastern  Hemisphere 
location  (34  degees  east  longitude). 

It  will  be  seen  that  the  transition  between  Phase  B and  Phase  C 
will  not  constitute  a major  change  in  the  type  of  activities  being 
conducted,  rather  only  the  level  and  extent  of  most  activities  will  be 
reordered. 

A typical  day  during  this  time  period  will  consist  of: 

Ground  tracking  systems  (Laser,  C~Band  and  Doppler)  data  collections 
activities  on  a global  basis  to  the  maximum  extent  possible  commensurate 
with  power  budget,  other  systems  activities,  and  in*'estlgator  needs. 

Altimeter  experiment  data  collection  activities  commensurate 
with  power  budget,  investigator  needs  and  globally  within  the  constraints 
of  the  ground  TM  station  coverage. 
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In  addition  to  this  major  activity,  intermittent  calibration  and 
evaluation  activities  will  be  conducted  at  a nominal  rate  of  once 
per  month, 

PHASE  D - Global  Activities 

Phase  P activities  will  be  conducted  during  the  time  in  which  ATS-P 
is  being  maneuvered  from  the  Western  to  the  Eastern  Hemisphere. 

During  this  period,  four  arcs  of  SSE  data  per  day  will  be  scheduled 
for  a total  of  160  5S£  data  arcs*  Each  arc  will  be  about  45  minutes 
long  (i.e.,  compatible  with  the  period  of  mutual  visibility  between 
ATS-F  and  GEOS-C). 

The  SSE  arcs  will  be  scheduled  to  provide  a complete  5^  X 5®  grid 
pattern  within  the  ATS-E  coverage  areas.  This  will  be  accomplished  by 
selecting  two  consecutive  ascending  and  two  consecutive  descending  passes 
par  day  (spaced  about  12  hours  apart  so  that  they  cross  near  the  equator), 
and  basically  following  these  passes  each  day  until  the  grid  is  complete. 
Arc  selection  as  described  above,  will  effectively  accomplish  the  5 degree 
gridworlc  since  GEOS-C  ground  traces  are  offset  by  about  5 degrees  to  the 
eaf.t  per  day. 

In  addition.  Interface  with  ATS-F  during  the  entire  40-day  period 
should  allow  a redundant  5- degree  grid  pattern  (offset  by  about  2 degrees 
to  the  east)  to  be  placed  In  an  area  extending  about  100  degrees  in  long- 
itude, This  additional  5-degree  pattern  will  be  placed  approximately 
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midway  within  the  ATS-F  coverage  area  (e..e,  » between  approximately 
80  degrees  west  longitude  and  20  degrees  east  longitude) • 

During  each  of  the  SSE  arcs  described  above,  radar  altimeter  data  will 
be  scheduled  over  most  portions  of  the  arcs  which  are  over  ocean  areas. 
Therefore,  a 5-degree  gridwork  of  altimetry  data  will  be  collected 
simultaneously  with  the  SSE  data*  In  addition,  C-Band,  laser,  and 
doppler  data  will  be  scheduled  to  support  orbit  determination  for  the 
SSE  and  radar  altimeter  data  arcs  scheduled  during  this  period. 

PHASE  E - Localized  Grid  Densification 

Phase  E activities  will  begin  at  the  time  when  activities  associated 
with  the  ATS-F  have  been  completed.  It  is  expected  that  Phase  E 
activities  will  be  essentially  the  same  as  those  described  in  Phase  C 
except  that  during  the  latter  time  period  the  S-Band  ground  tracking 
network  stations  should  all  be  modified  compatible  with  the  GEOS-C 
instrumentation  and  will  assume  a more  active  role. 


INVESTIGATION  PLANS 


Introduction 

In  support  of  the  G20S-C  Mission  objective.-^  proposals  were  solicited 
in  thirteen  specific  investigation  categories  with  provision  for  investi- 
gations in  other  categories  if  these  could  be  shown  to  be  compatible  with 
long  range  Earth  and  Ocean  Physics  Application  Program  objectives.  The  thirteen 
specific  investigation  categories  can  be  summarized  as  follows: 

Ocean  Gooid  Determination 

This  proposal  category  includes  all  proposals  for  the  determination 
of  the  geometry  of  mean  sea  lev.l  using  altimetry  data  alone  or  in  combin- 
ation with  other  data  types.  The  satellite  altimeter  observations  will 
provide  measurements  of  the  height  of  the  satellite  above  the  ocean  surface. 

This  data  can  be  used  directly  to  estimate  the  ocean  geold,  provided  the 
satellite  position  can  be  determined  with  sufficient  accuracy  and/or 
errors  in  satellite  position  corrected  for.  Investigations  in  this 
category  may  call  for  the  combination  of  altimeter  information  with  geo id 
information  obtained  from  existing  surface  gravimetry,  satellite  gravity 
field  information,  deflection  of  the  vertical  information  and  geocentric 
station  position.  The  computation  of  improved  satellite  gravity  fields 
and  station  positions  should  not  be  included  within  this  investigation 
category. 

One  of  the  Important  results  expected  to  be  obtained  from  the  GEOS-C 
altimeter  is  improved  definition  of  the  ocean  geoid.  At  present  worldwide 
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knowledge  of  the  ocean  geoid  is  only  available  from  satellite  gravity  field 
data  which,  at  best,  define  variations  with  widths  of  the  order  of  1500 
kins  or  larger.  Over  restricted  areas  of  the  world,  where  dense  surface 
gravity  data  is  available,  detailed  geoids  can  be  computed  defining  geoid 
variations  with  widths  down  to  100  kms  or  smaller.  These  detailed  geoids 
in  local  areas  have  demonstrated  that  geoid  variations  of  10  to  20  meters 
are  commonly  generated  by  the  wavelengths  of  less  than  1500  kms  and  are 
not  present  in  the  satellite  fields*  It  has  also  been  shown  that  even 
wavelengths  of  less  than  100  kms  can,  at  times,  produce  variations  of  up 
to  10  meters.  Therefore,  the  satellite  altimeter  with  precision  and/or 
accuracy  of  1 to  2 meters  has  the  potential  for  greatly  Increasing  our 
knowledge  of  the  ocean  geoid  in  those  substantial  parts  of  the  ocean 
where  no  detailed  surface  gravity  erists  as  well  as  contributing  to  in- 
creased accuracy  in  those  areas  where  surface  gravity  and  other  types  of 
gravity  data  exist. 

Determination  of  geoid  heights  from  altimeter  measurements  requires 
that  the  altimeter  measurements  be  reduced  in  conjunction  with  satellite 
orbit  information.  Questions  exist  as  to:  (1)  the  uesL  method  of  redact ion 
to  eliminate  possible  systematic  orbit  and  altimeter  errors,  (2)  the  best 
set  of  parameters  to  represent  the  geoid,  and  (3)  the  best  method  for 
combination  of  altimeter  data  with  other  types  of  data  for  improved  geoid 
determination.  Ten  GEOS-C  investigations  were  proposed  in  the  geoid 
computation  area.  These  Investigations  represt'ut  a number  of  approaches 
to  geoid  determination  from  altimeter  data  and  are  expected  to  provide 
answers  to  the  questions  posed  above.  In  order  to  answer  the  questions 
posed,  evaluations  must  be  made  of  the  actual  altimeter  results.  Several 
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investigations  have  as  their  objective  the  carrying  out  of  comparisons 
of  results  with  external  standards  in  order  to  make  such  evaluations. 

Ocean  Tides 

Repeated  measurements  ove^;  a section  of  the  ocean  using  the  GEOS->C 
altimeter  has  the  potential  for  allowing  determination  of  the  time 
variable  effects  of  the  tidal  attractions  of  the  sun  and  moon  on  sea 
surface  topography.  At  present  most  measurements  of  ocean  tides  are 
made  at  coastal  stations  where  the  tidal  effects  are  strongly  Influenced 
by  local  bathymetric  effects.  Although  several  theories  exist  which 
permit  theoretical  computation  of  deep  ocean  tides  only  limited  numbers 
of  measurements  of  deep  ocean  tides  have  been  made,  utilizing  bottom  tide 
meters. 

The  satellite  altimeter  has  the  potential  for  rapid  global  deter- 
mination of  ocean  tides.  To  be  of  maximum  use,  accuracies  of  the  order  of  10cm 
will  be  required  for  altimeter  derived  tidal  measurements.  However,  GEOS-C 
even  with  the  lesser  accuracy  of  its  altimeter  should  allow  evaluation 
of  various  techniques  for  recovery  of  tide  data  from  satellite  altimeter 
measurements.  To  be  of  maximum  valu^  tidal  analyses  of  GEOS-C  altimeter 
data  must  be  carried  out  In  areas  where  ground  truth  In  the  form  of 
bottom  tide  meter  data  is  available.  Five  Investigations  In  this  category 
were  proposed. 

Sea  State  Determinations  / 

In  addition  to  giving  the  distance  between  the  spacecraft  and  the 
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ocean  surface » the  GEOS-C  altimeter  datai  through  analysis  of  the  character- 
istics of  the  return  pulse,  is  expected  to  provide  Information  on  the 
sea  state.  In  particular,  Information  on  mean  wave  height,  wave  period, 
and  wave  propagation  direction  may  be  determinable.  Although  theoretical 
studies  and  aircraft  radar  altimeter  data  analyses  have  been  carried  out, 
considerable  effort  is  needed  to  determine  the  degree  to  which  various 
types  of  sea  state  data  can  be  extract  -4  from  a satellite  altimeter  and 
to  identify  the  best  methods  for  carrying  out  extraction  of  tha  Information, 
The  bulk  of  the  ten  investigations  proposed  on  GEOS-C  for  sea  state 
determination  analyses  are  aimed  at  evaluation  of  feasibility  and  iden- 
tification of  best  methods  through  comparis'^n  of  results  obtained  from  the 
GEOS-C  altimeter  with  ground  truth  Information  on  sea  state  and  with  data 
obtained  from  aircraft-borne  radar  instruments.  In  addition  to  analysis 
of  GEOS-C  data  in  terms  of  sea  state  parameters,  the  ohjectlvos  of 
investigations  include  development  of  Information  for  use  in  the  design 
of  future  satellite  radar  altimeters  and  determli  ation  of  potential  bias 
introduced  Into  altimeter  sea  surface  topograpny  determinations  due  to 
sea  state. 

Quasl-Statlonary  Departures  from  the  Marine  Geold 

This  proposal  category  includes  all  altimete-  analyses  designed  to 
Investigate  non-periodic  deviations  of  sea  level  from  an  equlpotentlal  sur- 
face* It  also  includes  analyses  of  altimeter  data  to  determine  sea 
slopes  associated  with  such  phenomena  as  currents  and  wind  setup.  It 
does  not  Include  analyses  relating  to  wave  pehnomena,  ocean  tides,  or 
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investigations  directed  specifically  to  the  determination  of  the  geoid. 

The  sea  surface  topography  which  will  be  measured  by  the  GEOS-C 
altimeter  is  a function  primarily  of  variation  of  the  fo  ce  of  gravity 
over  the  earth's  surface » changes  in  atmospheric  pressure  from  point  to 
point  on  the  ocean  surface  surface,  density  structure  of  the  water 
column,  surface  wind  effects,  dynamical  effects  due  to  ocean  currents, 
and  tidal  effects.  If  only  gravitational  forces  (Including  rotation) 
were  present,  the  sea  surface  topography  would  be  coincident  with  the 
geoid.  The  effects  of  atmospheric  pressure  variations,  ind  forces 
and  tides  are  time  variable  with  a reasonably  high  temporal  frequency. 

The  effects  of  density  structure  of  the  water  column  and  currents  are  usually 
considered  to  be  quasJ -stationary  departures  from  the  ..  ^ ',  even  though 
the  effects  of  currents  do  shift  over  restricted  areas  :he  surface. 

One  of  the  primary  aims  of  the  Earth  and  Ocean  Physics  Applications 
Program  Is  to  determine,  from  altimeter  measurements,  departures  of 
sea  surface  topography  from  the  marine  g4old  Oms  to  wa^'-r  motion.  The 
reason  for  this  li.terest  lies  in  the  fact  that  the  velocity  and  volume 
of  water  in  motion  can  be  infered  from  these  departures.  Although  the 
GEOS-C  altimeter  is  not  expected  to  be  accurate  enough  to  provide  in- 
formation which  is  scientifically  useful,  it  is  expected  to  be  accurate 
enough  to  test  out  the  concept  and  permit  development  of  methodologit^s 
which  can  be  used  with  data  from  latter,  more  accurate,  satellite  altimeters. 
Two  Investigations  are  planned  in  this  category. 

Gravity  Model  Improvement 


This  proposal  category  includes  all  analyses  of  GEOS-C  altimeter  and 


tracking  data  whose  ultimate  objective  is  the  determination  of  an  improved 
earth  gravity  field.  These  include  both  normal  perturbation  analyses 
combining  GEOS-C  tracking  data  from  other  satellites  and  analyses  in 
which  the  altimeter  geold  height  information,  SSE  rate  information  or 

other  tracking  data  are  combined  with  existing  Information  for  gravity  field 
improvement. 

Improvement  of  the  existing  gravity  models  is  required  to  achieve 
EOPAP  Program  goals  from  three  viewpoints.  First,  satisfaction  of  a 
number  of  EOPAP  goals  requires  improved  satellite  orbit  determination 
which  is,  to  a large  extent,  dependent  upon  an  improved  gravity  field. 
Second,  determination  of  effects  of  ocean  currents  on  sea  surface 
topography  requires  high  accuracy  geoidv  with  which  altimeter  derived 
sea  surface  topography  can  be  compared.  Increased  geold  accuracy  requires 
Increased  accuracy  in  knowledge  of  the  gravity  field.  Finally,  inter- 
pretation of  an  improved  gravity  field  offers  the  potential  of  Increased 
understanding  of  platetectonics  and,  therefore,  of  the  mechanisms  producing 
earthquakes. 

Gravity  field  information  can  be  derived  from  GEOS-C  in  three  ways: 

(1)  By  comDlning  information  on  the  perturbations  of  GEOS^C  from  tracking 
data  with  data  from  other  satellites  in  a general  perturbation  analysis; 

(2)  By  analysis  ^of  satellite- to-sat ellite  tracking  data  in  the  same  manner 
as  Lunar  Orblter  and  Apollo  data  was  analyzed  to  obtain  residual  llne-of- 
Bight  accelerations  or  compatible  gravity  anomaly  information.  Siv  GEOS-C 
investigations  fall  within  these  areas.  Their  proposed  Investigations  often 
include  the  combination  of  GEOS-C  data  with  gravity  field  information  from 
other  sources* 
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CeoloRical  Investigations 


One  important:  use  of  the  geoid  results  to  be  derived  from  the  GEOS<-C 
altlmetev  data  will  be  Interpretation  in  terms  of  the  geological  and 
geophysical  significance  of  the  results.  The  GEOS-C  altimeter  results 
can  be  of  particular  value  in  extending  Informatlo.i  to  areas  where  little 
or  no  surface  gravity  presently  exists.  Two  investigations  were  proposed 
in  this  category. 

Solid-Earth  Dynamics 

This  proposal  category  Includes  all  analyses  Involving  the  determination 
of  earth  tides,  polar  motion,  and  changes  in  rotation  rate  of  the  earth. 

It  also  includes  determination  of  very  precise  positions  on  the  earth's 
surface  using  GEOS-C  tracking  data  for  such  purposes  as  determination 
of  fault  notion  and  crustal  plate  motion. 

High  precision  tracking  of  the  GEOS-C  satellite,  particularly  by  the 
submeter  precision  laser  sy terns,  will  allow  derivation  of  improved  infor- 
mation on  the  dynamics  of  the  solid  earth.  Determinations  can  be  made  of 
the  gravitational  and  geometric  effects  of  solid  earth  tides  and  of  the 
motions  of  the  earth's  pole  Including  chandler  motion,  yearly  motion,  and 
the  diurnal  wobble.  Four  invesitgations  were  proposed  in  this  category. 

Interconparison,  Evaluation  and  Calibration  of  Instrumentation  Systems 

This  investigation  category  Includes  all  investigations  whose 
objective  is  the  evaluation  and  calibration  of  altimeter,  sate 111 te-to- 
satelllte  tracking  and  ground  tracking  instrumentation  to  be  used  with 
the  GEOS-C  mission.  This  includes  both  evaluation  of  the  on-board 
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instrumentation  and  the  ground  systems.  In  this  investigation  category 
are  placed  all  instrument  intercomparison  investigations  and  studies 
related  to  instrumentation  technology. 

The  GEOS-C  Project  will  undertake  the  primary  evaluation  and 
calibration  activity  with  respect  to  the  radar  altimeter.  Six  investigator 
are  planned  to  be  associated  with  evaluation  and  calibration  activities 
relative  to  ground  tracking  instruments  and  the  satellite-to-satellite 
experiment.  Calibration  and  evaluation  of  ground  tracking  instruments 
can  be  carried  out  by  ^ ntercomparisoa  analyses.  These  can  Involve  both 
co-located  inter comparisons  and  inter comparisons  involving  reference  orbits. 
Since  the  satellite-to-satellite  experiment  Involves  new  Instrumentation, 
special  emphasis  will  be  given  to  evaluation  and  calibration  of  these 
results. 

Ground  Truth  termination 

This  Investigation  category  Includes  all  investigations  whose  objective 
is  the  collection  of  data  from  ground,  ship,  and  aircraft  based  systems 
and  the  use  of  this  data -to  evaluate  the  charsct«r-t  cs  of  the  satellite 
systems. 

In  order  to  calibrate,  evaluate,  and  utilize  data  taken  by 
instruments  on  board  the  GEOS^C  spacecraft,  it  will  be  necessary  to  have 

available  certain  types  of  ground  truth  for  comparative  purposes.  Three 
Investigators  have  proposed  to  provide  some  of  this  ground  truth  data. 
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Tracklag  Station  Location  Improvement 
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This  Investigation  category  includes  all  investigations  vhose  primary 
objective  is  the  determination  of  the  location  of  tracking  stations 
where  the  objective  is  geodetic  in  nature  and  is  not  for  earth  dynamic 

r 

purposes. 

A number  of  types  of  tracking  data  taken  using  the  GEOS-C  satellite 
can  be  used  to  provide  improved  station  location  information  which  will  I, 

be  useful  in  support  of  altimeter  calibration  and  to  support  other  project 
objectives*  GEOS-C  will  provide  data  from  new  stations,  data  of  higher  - 

accuracy  than  previously  available,  and  data  from  new  instrumentation  ; 

types  such  as  VLBl  measurements.  Five  investigations  have  been  proposed 
in  this  category.  j 


Orbit  Determination  Improvement 

This  proposal  category  includes  all  investigations  whose  end  objective 
Is  orbit  determination  Improvement.  Indirectly,  GEOS-C  can  be  expected 
to  support  Improved  orbit  determination  by  providing  improved  gravity 
field  information.  However,  this  category  will  emphasize  new  types  of 
tracking  information  such  as  the  SSE  and  altimeter  data  and  Its  capability 
to  support  improved  orbit  determination.  Two  investigators  propose  to 
study  the  use  of  GEOS-C  data  for  this  purpose. 

Data  Management /Informat ion  Processing 

This  proposal  category  includes  investigations  whose  objective  is 
the  development  of  methods  and  techniques  for  managing  and  processing  the 
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data  taken  by  the  various  ins t rumen f^t ion  on  the  GEOS-C  spacecraft. 

This  includes  the  development  of  data  editing  and  pre-processing  techniques. 
Specifically,  investigations  are  directed  toward  those  systems  expected 
to  be  most  useful  in  future  earth  and  ocean  physics  applications  activities 
and  involve  advanced  techniques  applicable  to  future  activities. 

One  investigation  falls  in  the  general  category  rf  data  management 
and/or  information  processing  relative  to  the  altimeter. 

Unique  System  Investifijatlons 

Certain  GEOS-C  investigations  are  proposed  which  are  uniquely 
associated  with  a particular  instrumentation  and  do  not  fit  into  any  of 
the  preceeding  twelve  categories.  Two  Investigations  are  proposed  in 
this  category  and  deal  with  the  C-Band  and  Altimeter  systems. 
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CALIBRATION  PASSES  FOR  ALTIMETER  STABILITY  EVALUATION 


C-BAND  RA33AR  PARIICIPATION  REQUIKEi 


rESAZ  ONES  Azores  1-6  months 
'S-36  DFVLR  Germany  1-6  months 
'S-16  RAE  England  1—6  months 
S - 36  WSMR  Green  River,  Utah  1 - 6' months 


GEOS-C 

altimeter  system  characteristics 


CHARACTERISTIC 

LONG  PULSE  MODE 

S110T:T  PULSE  MODE 

OUTPUT  TUBE 

MAGNETRON 

TRAVELING  WAVE 

OUTPUT  POWER 

2 Kv 

2.5  K.W 

PULSE  WIDTH 

200  N SEC. 

12  N SEC, 

PRF 

100  PULSE  BURSTS/SEC. 
OF  16  PULSES/BURST 

100  PPS 

ACQUISITION  TIME 

<5  SEC. 

<5  SEC. 

OUTPUT  (ALTITUDE) 

32  BIT  AVERAGE 

32  BIT  AVERAGE 

NOMINAL  AVERAGE  TIME 

1 SEC. 

0.1  SEC, 

RANGE  RESOLUTION 

6.25  N SEC. 

1,56  N SEC. 

RANDOM  ERROR 
CORRELATION 

<0.05 

<0.33 

PRECISION 

<■1  m 

<0.6  m 
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APPENDIX  L 


GEOS-C  C-BAND  OPERATIONS/SUPPORT 


By 

BEN  JACKSON 

NASA/Wallops  Station 
Wallops  Island,  Virginia  23337 
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<;EaS-C  C-ftANP  SCHEDULING  S DATA  HANDLING  PROCEDURES 
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STAtJIEY/OPERATE  MODE  SWITCHING. 
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■ DEFENSE  MAPPING  AGENCY  TEST 

OBJECTIVES  FOR  GEOS-C 


By 

MAJ.  LARRY  BEERS 

Defense  Mapping  Agency 
Bldg,  56,  US  Naval  Observatory 
Washington,  DC  20305 


m 


SST  DATA  DURING  ATS  DRIFT  WOULD  NOT  PROVIDE  DENSm  REQUIRED  AND  mJlL  BE 
ACQUIRED  TOO  LATE  IN  THE  MISSION  TO  SATISFY  THIS  REQUIREMENT 
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APPENDIX  0 


GEOS-C  C-BAND  WORKING  GROUP 


7 MARCH  1974  MEETING 


Vandenberg  AFB,  California  93437 


AGENDA 


GEQS-C  C-Band  Working  Group 
March  7,  1974  - 1330^630  PPT 


A.  C-Band  Experiment  Sutiwary 

B.  Sunmatlon,  C-Band  Systems  Objectives 

C.  Pre-Launch  Objectives 

D.  Quick-Look  Activities 

E.  Normal  Post-Launch  Objectives 

F.  Scheduling  and  D^ta  Handling  Procedures 

6,  Coherent  C-Band  Transponder  (Vega  Model  355C) 
H.  Discussion 
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6E0S-C  C'Band  Working  Group 


Name 

Agency/Organization 

Address 

Telephone 

ARLOWE,  H.  Duane 

Sandia  Labs. 

Div.  9474,  P.O.  Box  5800 
Albuquei'que,  NM  87115 

505-264-6556 

^EERS,  Larry  D.,  Haj. 

HQ  OMA/PRA 

Bldg  56,  Naval  Observatory 
Washington,  DC  20305 

202-254-4455 

>6£LGIN,  John  A. 

PMR 

Range  Operations,  Code  3274 
Point  Mugu,  CA  93041 

805-982-8139 

BERBERT»  John 

NASA/GSFC 

Code  932 

Greenbelt,  HD  20771 

301-982-5055 

♦BORREGO,  Arturo 

USHR 

STEWS-IO-E 

White  Sands  Missile  Range, 
NM  88002 

915-678-3220 

CARNEY,  DeVere 

RCA  Service  Co. 

8855  Annapolis  Road 
Lanham,  MO  20601 

301-552-9600 

CHIN,  Ball 

WSMR 

STEWS-ID-E 

White  Sands  Missile  Range, 
NM  88002 

915-678-5405 

DEMPSEY,  Don 

RCA/MSRD 

Bldg  lOP-239 
Moorestown,  NJ  08057 

609-963-8000/ 

609-S63-3283 

EBAUGH,  Larry 

AFETR  HQ 

Attn:  DOOT 

Patrick  AFB.  FL  32925  i 

305-853-5941/ 

AV467-5941 

GOBLE,  Phil 

' SAHTEC/FEL 

P.O.  Box  1886,  Bldg  8310 
Vandenberg  AFB,  CA  93437 

805-276-4189 

GRAVES,  Kenneth 

SAMTEC/ROOE 

Bldg  7000 

Vandenberg  AFB,  CA  93437 

AV276-8995 

GREENE.  John  A.,  MaJ. 

AFETR/ RCA 

Bldg  389 

Patrick  AFB,  FL  32925 

305-494-2191 

HALL,  Edgar 

SAMTEC/FEC 

P.O.  Box  1886,  Bldg  8510 
Vandenberg  AFB,  CA  93437 

805-866-7578 

AV276-7578 

HAWKINS.  William  E. 

j 

KASA/G3FC 

Code  861.3 
Greenbelt,  MO  20771 

301-962-6856 

♦HILLHOUSE,  Hilton 

AFETR  HQ/DOS 

Patrick  AFB.  FL  32925 
1-3 

305-494-7735 

Agency/Organi zati on 


SAFSCOM/KMRD 


KENNEDY,  Ootin  M. 

KRABILl,  William 
LANAN,  K.  F.,  L/Col. 
LUBAR,  Bertram  H. 
*MANNING,  Walter 

MITCHELL,  R.  D. 

MORGAN,  Donald  G. 

OSTlIE,  Dick 

PARKER,  Horace 

PATRICK.  John  \ 

*PAULUS,  William  K. 

PEASE,  C.  E. 

*RAD0M,  Stanley  R. 
RELF,  Kenneth  E.,  Dr. 
ROLLINS,  Clarence  M. 

ROY,  N.  A. 


SAHTEC/XPP 

AFETR/ENLT 

AFETR/RML 

RCA/MSRD 

Vega  Prec.  Labs 

SAHTEC/FEE 

RCA  Service  Co. 

Sandia  Labs. 


SAHTEC/ROPP 
SAHTEC  HQ/CA 
AFETR/RCA/RWL 
AFETR/ENIL 


Cmdr,  Safeguard  System  Command 

205-895-453J 

SSC-RO/T.  Keeney 
P.O.  Box  1500 
Huntsville,  AL  35B07 

Patrick  AFB,  FL  32925 

305-494-5129/ 

AV854-5034 

Wallops  Island,  VA  23337 

804-824-3411 

Vandenberg  AFB,  CA  93437 

805-866-7875 

Patrick  AFB,  FL  32925 

AV659-7096 

Bldg  981 

Patrick  AFB,  FL  32925 

305-494-2211 

Bldg  108-239 

609-963-5000/ 

Moorestown,  NO  06057 

609-963-3715 

600  Foil In  Lane 
Vienna,  VA  22180 

703-938-6300 

P.O.  Box  1886,  Bldg  8310 
Vandenberg  AFB,  CA  93437 

j AV276-3436 

8855  Annapolis  Road 
Lanham,  20801 

301-552-9600 

Tonopah  Test  Range 
P.O.  Box  871 
Tonopah,  NV  89049 

702-985-6277 

Div.  9426 

Albuquerque,  NM  87115 

505-264-1664 

Vandenberg  AFB.  CA  93437 

AV276-7631 

Vandenberg  AF8,  CA  93437 

j AV276-7022 

Patrick  AFB.  FL  32925 

* 305-494-2950 

Attn:  ENIL 
Patrick  AFB,  FL  32925 

AVb59-4271 

Pocomoke  City,  MD  21851 

j 301-957-2484 

NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 
WALLOPS  STATION 
Wallops  Island.  Virginia  2333? 


REPLY  TO 
hi  Hi  OP; 


PMS  (GEOS-C) 


tm  2 1974. 


TO:  Distribution 

FROi^:  PMS/GEOS-C  C-Band  Subsystems  Manager 

SUBJECT:  GEOS-C  C-Band  Working  Group 


The  6E0S-C  C-Band  Working  Group  Meeting  was  '.eld  on  March  7,  1974,  at 
Vandenberg  Air  Force  Base,  California,  in  Building  7000  Theatre,  as  a 
part  of  the  SAMTEC  Conference  on  Coherent  Radars  for  Range  Instrumenta 
tion.  Since  most  of  the  cognizant  Department  of  Defense  C-Band  instru 
mentation  personnel  v/era  in  attendance  for  the  conference,  which  began 
on  March  5,  1974,  it  was  appropriate  to  conclude  the  conference  with 
the  C-Band  Working  Group  Meeting. 

It  was  easily  apparent,  during  the  entire  conference,  that  the  6E0S-C 
C-Band  instrumentation  would  be  of  invaluable  assistance  to  all  of  OOD 
as  a means  of  calibrating  and  evaluating  the  performance  and  accuracy 
of  their  tracking  systems.  Several  requests  were  maae  during  tne 
conference,  by  persons  not  originally  members  of  the  Working  Group, 
to  attend  this  session,  and  38  persons  (see  Enclosure  1)  were  present. 
Sandia  Laboratories  and  the  Defense  Mapping  Agency  v/ere  the  new 
agencies  represented. 

A general  sunr/nation  of  the  C-Band  mission  objectives  was  given  and 
general  dialogue  regarding  the  impart  of  5AMTEC  as  lead  DOO  range  on 
scheduling,  coordination,  and  data  handling  ensued.  The  basic  con- 
sensus of  opinion  by  NASA,  Wallops,  and  SAMTEC  was  that  NASA,  Wallops, 
would  do  all  advanced  sc.heduling  of  DOD  radars  through  SAITEC.  Any 
changes  to  the  schedule  would  be  done  directly  between  Wallops  and 
the  specific  range  involved  with  an  information  copy  to  SAMTEC.  The 
reverse  would  hold  true  for  the  schedule  items  originally  requested 
by  DOD  for  tracking  support  by  various  support  ranges  if  conflicts 
arose.  All  ranges  tracking  would  submit  a tracking  report  at  approxi- 
mately T+1  hour  via  TWX  to  Wallops  and,  if  the  tracking  was  in  support 
of  NASA,  Wallops,  use  the  mails  to  send  the  data  packet  directly  to 
'Wallops.  Questions  arose  concerning  distant  or  remote  sites,  and  it 
was  concluded  that  the  same  procedures  would  be  utilized  for  them. 

Questions  began  to  arise  coi. earning  the  lack  of  formal  information  at 
the  respective  support  ranges  for  the  6E0S-C  mission,  and  copies  of 
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the  Program  Introduction  (PI)  Document;  was  distributed  to  a representa- 
tive of  each  range.  This  document  formalizes  the  program  and  it  becomes 
an  official  support  requirement  for  each  range.  Suggestions  v/cre  made 
that  perhaps  a GEOS-C  team  be  invited  to  offer  presentations  similar  to 
those  in  the  SAMTEC  Agenda  to  the  Conmanding  Officer  at  each  range,  or 
that  an  invitation  be  extended  to  the  GEOS-C  team  to  make  a presenta- 
tion to  the  Range  Conmanders  Conference  in  the  near  future. 

Despite  the  factor  that  expenditures  for  direct  support  are  required 
by  submission  of  the  PI  according  to  a COD  directive  being  implemented, 
commencing  July  1,  1974,  most  persons  felt  that  the  DOD  requirements  on 
the  GEOS-C  mission  would  offset  any  costs  incurred  by  NASA. 

Many  topics  were  openly  discussed  and  explanations  offered  by 
H.  R.  Stanley,  W.  B.  Krabill,  and  A.  R.  Selser  of  NASA,  Wallops; 

W.  E.  Hawkins  and  J,  Berbert  of  NASA,  GSFC;  and  the  undersigned  on 
the  multiple  facets  of  the  GEOS-C  Program.  Listed  below  are  these 
topics  and/or  the  requirements  as  stated  by  the  range  representatives 
starred  on  Enclosure  1. 

a.  Sandia  Laboratories  at  the  Tonopah  Test  Range  requested  that 
they  be  permitted  to  con^nence  tracking  GEOS-B  in  preparation  for  the 
upcoming  GEOS-C  mission.  The  procedure  for  requesting  such  support 
was  explained,  as  well  as  the  available  power  remaining  on  GEOS-B, 
for  this  tracking.  Suggestions  ware  made  that  tiicy  coordinate  those 
tracks  for  acquisition  purposes  with  SAi'iTEC  and  WSMR. 

b.  The  ETR  must  solve  inhouse  the  interface  between  Range 
Measurement  Laboratory  (RI-IL)  and  the  remainder  of  the  range  such  that 
a single  scheduling  interface  be  made  through  L.  Ebaugh,  AFETR  Head- 
quarters. 


c.  Some  arrangement  should  be  made  by  NASA  to  identify  all 
support  requests  as;  REQUIRED,  DESIRED,  or  MISSION  CRITICAL.  It 
was  further  suggested  that  Major  L.  Beers/DiiA  attempt  to  establish 
military  priorities  for  C-Band  support. 

d.  The  required  NASA  pre-,  setup-,  and  post-calibrations  for 
data  collection  will  be  made  to  DOD  and  included  in  the  Operations 
Requirements  (OR)  Document  as  specifications  for  operations  in  support 
of  GEOS-C,  This  is  essential  for  ranges  having  multiplicities  of  the 
same  type  of  radars,  since  each  radar  has  different  characteristics. 
Raw  range  data  should  be  forwarded  to  Wallops. 

e.  The  proposed  14-day  advanced  schedule  requirements  were 
requested  to  be  made  30  days  to  allow  advance  planning  by  the  support 
ranges. 
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f.  ^ An  ranges  have  the  capability  of  generating  their  own  angles, 
if  provided  either  the  NORAD  or  BROUWER  elements,  BROUWER  elements  will 
be  provided  by  the  Project  through  GSFC. 

g.  A series  of  coordinated  tracks  will  be  set  up  utilizing  GEOS-B 
to  simulate  a GE03-C  type  operation  including  data  handling.  This  may 
have  to  te  done  in  segments  in  view  of  the  available  power  for  GEOS-B. 

Ii.  Many  of  the  ranges  have  altered  or  added  to  their  basic  requlre- 
ments/objeotives  for  the  GEOS-C  mission  and  are  to  supply  the  undersigned 
these  updated  requirements/objectives  by  the  end  of  March,  A short 
summation  of  each  is  indicated  below. 

(1)  ETR  - C-Band  data  for  self-calibration  and  coordination 
with  SAf-lJEC  for  USNS  ARNOLD.  Laser,  and  possibly  Doppler-  data  is 
requested  from  collocated  sites, 

(2)  WSMR  - C-Band  data  for  calibration  standards  and  inter- 
comparisons with  other  ranges. 

(3)  PMR  " C-Band,  Doppler,  and  Laser  data  coordinating  with 
other  ranges,  (Doppler  needs  were  passed  on  to  DMA  for  resolution.) 

(4)  W1R  - C-Band  data  for  self-calibration,  Doppler  orbits, 

TM,  and  altimeter  data  are  requested.  The  ALCOR  radar  was  mentioned, 
but  efforts  are  being  made  by  NASA  not  to  utilize  this  system,  (TM, 
Doppler,  and  Altimeter  needs  were  passed  on  to  OfiA  for  resolution.) 

(5)  Sandia  - C-Band  data  for  improved  station  location  and 
calibrations  with  other  ranges  (SAIiTEC,  WSMR). 

(6)  SAMTEC  - C-Band,  Altimeter,  and  Laser  data  for  calibra- 
tion coordination  with  other  ranges. 

Much  interest  has  been  generated  within  the  group,  and  the  unanimous 
consensus  of  opinion  is  that  these  Working  Group  Meetings  should  be 
held  every  two  months  throughout  the  lifetime  of  the  GEOS-C  mission, 
possibly  on  a range-host  rotacable  basis. 

Earl  B.  Jacl^Son 
' Enclosure 
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